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I.    INTRODUCTION 

A.        BACKGROUND 

As  costs  for  material  and  energy  continue  to  rise,  the 
need  for  smaller,  more  efficient  heat  exchangers  continues 
to  grow.  For  those  applications  involving  marine  vehicles, 
size  and  weight  limitations  dictate  the  use  of  small, 
highly-efficient  heat  exchangers.  Reducing  the  size  and 
weight  of  steam  condensers  used  aboard  U.  S.  Navy  ships  for 
propulsion  and  electrical  power  generation  would  result  in 
lower  material  costs,  and  help  to  aleviate  the  cramped 
conditions  so  typical  of  machinery  spaces.  Thus,  there  is  a 
strong    motivation   for  continued  study  in    this   area. 

The  effectiveness  of  condensers  is  limited  by  the  water- 
side, vapor-side,  and  wall  thermal  resistances  of  the 
condenser  tubes.  Reducing  any  one  of  these  will  contribute 
to  increased  hea t- transfer  performance,  and  smaller  physical 
size  of  condensers.  Methods  of  enhancing  the  vapor-side 
coefficient  include  the  use  of  "roped"  tubes,  fluted  tubes, 
drainage  strips  attached  on  the  tubes,  finned  tubes,  and 
coatings   applied      to    enhance      dropwise  condensation.  This 

thesis    concentrates    soley   on  finned    tubes. 

Since  the  late  1940s,  externally-finned  tubes  have  been 
used  to  increase  the  vapor-side  heat-transfer  coefficient  of 
tubes  used  in  refrigeration  systems;  but  condensers  used  in 
steam  systems,  such  as  shipboard  propulsion  plants,  continue 
to  use  smooth  tubes.  The  high  surface  tension  of  water, 
which  leads  to  its  tendency  to  flood  the  area  between  fins, 
has  resulted  in  a  widely-held  belief  that  such  tubes  are 
inappropriate  for  use  in  steam  systems.  Recent  studies 
[1,2],  however,  have  shown  that  finned  tubes  can  signifi- 
cantly  enhance    the    heat-transfer    rates   in    such    systems. 
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The  theoretical  treatment  of  condensation  on  finned 
tubes  is  extremely  complex  due  to  the  large  number  of  vari- 
ables and  physical  mechanisms  involved.  The  interaction  of 
gravitational  and  surface-tension  forces  lead  to  complex 
three-dimensional  flow  patterns,  which  are  further  dependent 
on  fin  spacing,  height,  and  thickness.  Other  variables 
include  heat  flux,  vapor  shear,  tube  diameter,  fin  shape  and 
fluid  properties  just  to  name  a  few.  In  view  of  the  above, 
any  theoretical  models  will  require  numerous  simplifying 
assumptions,  and  require  complex  computer  solutions 
involving  implicit  numerical  techniques  [3].  To  confirm  the 
validity  of  theoretical  models,  reliable  experimental  data 
which  cover  a  wide  range  of  relevant  parameters  must  be 
obtained.  '"he  availability  of  such  data  may  then  lead  to 
the  development  of  fairly  simple  experimental  correlations, 
which  could  be  used  in  the  design  stage  to  predict  the  heat- 
transfer  performance  of  finned  tubes. 

This  thesis  effort  is  a  continuation  of  research  being 
conducted  at  the  Naval  Postgraduate  School  (UPS)  under  a 
grant  from  the  National  Science  Foundation.  The  basic  test 
apparatus  used  to  collect  experimental  data  was  ruilt  by 
Krohn  [4].  Graber  [5]  provided  the  majority  of  instrumenta- 
tion, and  took  preliminary  data  as  the  system  experienced 
problems  with  non-condensing  gases  and  partial  dropwise 
condensation  on  copper  tubes.  Poole  [6]  made  further 
improvements  on  the  apparatus  as  well  as  on  the  instrumenta- 
tion and,  most  importantly,  assured  a  leak-free  apparatus. 
Unfortunately,  he  did  not  have  sufficient  time  to  produce 
useful  data,  mostly  due  to  the  considerable  time  spent  in 
systematically  locating  and  fixing  leaks,  and  due  to  the 
partial  dropwise  condensation  problem  that  had  not  been 
solved.  Using  this  system,  Georgiadis  [7]  was  finally  able 
to  obtain  complete  filmwise  condensation  on  copper  tubes. 
He   took  data   on   a  number   of  finned   tubes   with  fins   of 
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rectangular  shape,  as  well  as  on  smooth  tubes.  The  repeat- 
ability of  data  obtained  by  Georgiadis  demonstrated  the 
accuracy  of  the  test  apparatus  and  associated  instrumenta- 
tion vhich  was  used  essentially  without  modification  for 
this  investigation. 

The  overall  objectives  of  the  present  program  at  NPS 
includes  the  testing  of:  (a)  tabes  with  rectangularly- 
shaped  fins  to  find  the  best  fin  spacing,  thickness,  and 
height,  (b)  tubes  with  various  fin  shapes  to  find  the  best 
geometry  that  will  maximize  heat  transfer,  (c)  tubes  with 
different  tube-metal  thermal  conductivity,  and  (d)  the 
effect  of  vapor  shear  on  finned  tubes.  Georgiadis  [7] 
tested  a  total  of  25  tubes,  which  were  divided  into  five  fin 
spacings  (0.5,  1.0,  1.5,  2.0  and  4.0  mm),  five  fin  thick- 
nesses (0.5,  0.75,  1.0,  1.25,  and  1.5  mm),  and  two  fin 
heights  (1.0  and  2-0  mm).  At  the  conclusion  of  Georgiadis' 
test  program,  eight  new  tubes  with  fin  heights  of  0. 5  and 
1.5  mm,  fin  spacings  of  1.0,  1.5,  2.0,  and  4.0  mm,  and  a 
fin  thickness  of  1.0  mm,  remained  to  be  tested  to  complete 
the  sequence  of  tubes  with  rectangularly- shaped  fins.  The 
testing  of  these  tubes  and  additional  tubes  toward  objec- 
tives (b)  ,  (c)  ,  and  (d)  were  the  primary  goals  of  this 
thesis  effort  as  listed  in  the  next  section. 

B.   OBJECTIVES 

The  main  objectives  of  this  thesis  are  as  follow: 

1.  Take  data  on  tubes  with  rectangularly-shaped  fins  of 
various  fin  spacings  and  fin  heights  to  augment 
previous  data  [7 ], 

2.  Take   data  on   tubes  with   fins   of  different   shapes 

(triangular,  trapezoidal,  parabolic,  etc.), 

3.  Take  data  on  commercially-available  finned  tubes, 
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4.  Take      data      on      the      "optimum,"      rectangularly-shaped 
finned   tube   at   different    vapor   velocities, 

5.  Take   data    on      a   stainless   steel    tube      with    "optimum," 
rectangularly-shaped    fins,    and 

6.  Develop    a      preliminary   correlation    based    on      data    for 
tubes    with   rectangularly-shaped    fins. 


19 


II.  THEORETICAL  TREATMENT  OF  CONDENSATION  ON  HORIZONTAL 

FINNED  TOBES 

A.   FILH  CONDENSATION 

When  filmwise  condensation  of  a  vapor  takes  place  on 
smooth  tubes,  a  thin  layer  of  condensate  is  formed  which 
thickens  with  increasing  distance  around  the  perimeter  of 
the  tube.  This  condensate  layer  creates  a  thermal  resis- 
tance, which  can  limit  the  heat- transfer  performance  of  the 
tube.  This  film  thickness  and  its  thermal  resistance  can  be 
reduced  ty  the  use  of  external,  radial  fins  whish,  in  addi- 
tion to  a  surface  area  increase,  promote  surface-tension 
effects. 

In  1984,  Yau  et  al.  [1]  measured  the  enhancement 
provided  by  copper  finned  tubes  over  smooth  tubes  for  film- 
wise  condensation  of  steam.  Similar  experiments  by 
Wanniarachchi  et  al.  [2]  also  in  1984  confirmed  that  the 
observed  enhancements  were  greater  than  could  be  explained 
by  the  increased  surface  area  alone.  This  additional 
enhancement  may  be  a  result  of  the  surface-tension  forces 
which  act  to  thin  the  condensate  film.  Figure  2.1  schemati- 
cally depicts  this  phenomenon. 

The  effect  of  surface  tension  on  pressure  at  the  inter- 
face between  a  liquid  and  vapor  is  inversely  proportional  to 
the  radius  of  curvature  of  the  interface.  If  a  surface  has 
two  radii  of  curvature  at  right  angles  (i.e.,  r  and  r  ),  it 
can  le  shown  that: 
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o   =    surface    tension    of    condensate,    and 
A  F   =    the   prtssure   difference. 


Tin 


Condensate 


Figure  2.1    Schematic  of  Condensate  Profile  at  Dpper 

Tube  Surface. 


For  the  case  of  a  finned  tube,  the  radius  of  curvature 
around  the  outer  fin  perimeter  is  very  large  compared  to  the 
radius  of  curvature  around  the  fin  cross-section  profile 
(Figure  2.1).  Referring  back  to  aquation  (2.1),  it  can  be 
seen  that  the  smaller  radius  of  curvature  term  will  domi- 
nate, so  surface-tension  effects  around  the  fin  perimeter 
may    he    neglected. 

Because  of  the  convex  shape  of  the  condensate  film  at 
point  A,  the  pressure  within  the  film  at  this  point  is 
greater  than  the  surrounding  vapor  pressure.  In  a  similar 
manner,  the  pressure  at  point  B  is  less  than  the  surrounding 
vapor    pressure    owing    to    the      concave    shape    of   the    condensate 
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film.  The  relatively  flat  shape  of  the  condensate  film  at 
point  C  leads  to  a  pressure  essentially  egual  to  the 
surrounding    vapor    pressure.       These    pressures   are    given    by: 


where 


p     =  P     +  ^_  (2.2) 

A  v       rA 


%   =   Pv  "  77  (2-3) 


c         v  (2.4) 


Pv  =  surrounding  vapor  pressure, 

PA»  p  B*  pc  =  liquid  pressure  at  points  A,  B,  and  C,  and 

r  A/  r  q,    rc  =  radius  of  curvature  of  the  condensate 
film  at  points  A,  B,  and  C. 

At  point  B,  the  radius  of  curvature  is  small,  so  the  pres- 
sure at  point  B  is  less  than  the  pressure  at  point  C. 
Further,  from  equations  (2.2)  and  (2.3),  the  pressure  at 
point  A  is  greater  than  at  point  B.  In  reality,  the  pres- 
sure gradient  within  the  condensate  film  varies  along  the 
height  of  the  fin  due  to  the  continuously  varying  radius  of 
curvature  from  the  fin  tip  down  to  the  base  [3];  however,  to 
simplify  the  treatment  of  condensate  flow  toward  the  fin 
base,  the  pressure  differences  between  points  A  and  B  and 
points  C  and  B  may  be  simply  written  as: 
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APAB,  APCB   =  pressure  difference  between  points  A  and  B, 
and  points  C  and  B. 

As  can  he  seen,  these  pressure  differences  are  positive, 
resulting  in  condensate  flow  toward  point  B.  Due  to  the 
relatively  large  mass  of  condensate  at  point  5,  gravita- 
tional forces  dominate  over  surface-tension  forces  causing  a 
flow  of  condensate  around  the  tube  perimeter  at  the  fin 
base.  In  this  manner,  a  condensate  run-off  channel  is 
foraed  at  point  B,  and  the  improved  drainage  of  condensate 
thins  the  film  between  fins  and  ou  the  fin  surface.  This 
thinning,  in  turn,  reduces  the  thermal  resistance  through 
the  condensate  film,  thus  producing  an  enhanceaent  in  addi- 
tion to  the  gain  in  surface  area. 

It  should  be  stressed  at  this  point  that  the  shape  of 
the  condensate  film  is  highly  dependent  on  the  fin  geometry. 
Clearly,  then,  the  fin  shape  must  be  taken  into  considera- 
tion in  order  to  maximize  the  beneficial  effects  of  surface 
tension. 

The  above  gains  may  be  partially  or  totally  offset, 
however,  by  the  tendency  of  condensate  to  flood  the  area 
between  fins,    especially  on  the   lower  part  of   the  tubes. 
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The  extent  to  which  the  film  floods  the  tube  is  defined  by 
the  condensate  retention  angle  (  ip  )  ,  as  shown  in  Figure  2.2. 
In  the  flooded  portion  of  the  tube,  the  relatively  thick 
condensate  film  increases  the  thermal  resistance,  thereby 
leading  to  a  degradation  in  the  heat-transfer  performance. 
The  condensate  retention  angle  is  highly  dependent  on  fin 
spacing,  so  very  closely-spaced  fins  may  lead  to  completely 
flooded  tubes. 

A  number  of  studies  have  been  conducted  to  investigate 
the  performance  of  finned  tubes,  and  these  are  presented  in 
the  next  section. 

B.   CONDENSATE  RETENTION 

In  1946,  the  first  measurements  of  condensate  retention 
were  made  by  Katz  et  al.  [8].  These  measurements  were  made 
under  static  conditions  (i.e.,  no  condensation  taking 
place),  using  water,  aniline,  acetone,  and  carbon  tetrachlo- 
ride on  a  number  of  tubes  with  different  fin  densities  (276 
to  984  f  ins/m)  ,  and  fin  heights  (1.2  to  5.7  mm).  It  was 
shown  that  as  much  as  100£  of  the  tube  surface  could  be 
flooded  by  retained  condensate,  depending  mainly  on  the 
ratio  of  surface  tension  to  liquid  density  and  on  the  fin 
spacing. 

Eight  years  later,  Gregorig  [9]  recognized  that  surface- 
tension  effects  could  play  an  important  role  on  a  vertical, 
fluted  surface  which  uses  a  minimum  radius  of  curvature  at 
the  flute  tip,  which  gradually  increases  toward  the  trough. 
As  discussed  in  section  A  above,  this  variable  radius  of 
curvature  results  in  a  condensate  film  pressure  that 
decreases  toward  the  trough.  In  this  manner,  the  condensate 
layer  on  the  convex  region  of  the  flute  is  thinned  signifi- 
cantly, thus  pushing  the  condensate  into  the  troughs  (Figure 
2.3)  where  gravitational  forces  result  in  enhanced  drainage. 


24 


External  Diameter  of  fins 


Root  Diameter  of  fins 


Tube  Wall 


Retained  Condensate 


Rentention  Angle  ( C  ) 


Figure    2. 2 


Schematic    of    Condensate    Retention 
on    Finned   Tubes. 
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Since  the  heat- transfer  coefficient  is  inversely  propor- 
tional to  the  thickness  of  the  condensate  layer,  heat- 
transfer  performance  at  the  trough  is  degraded,  but  to  a 
lesser  degree  than  the  enhancement  experienced  by  the  crest. 
The  resulting  average  heat-transfer  coefficient,  therefore, 
shows  a  significant  enhance men t. 

THIN  FILM  REGION 


Figure   2.3        Cross   Section   of   the   Gregorig    Surface 


In  1981,  Rudy  and  Webb  [10]  made  measurements  of  conden- 
sate retention  angles  on  finned  tubes  with  three  different 
fin  densities  (748,  1024,  1378  fins/neter)  using  water, 
P.- 1  1  ,  and  n-pentane  as  the  working  fluid.  Later,  in  1983, 
Pudy  and  \lebh  [11]  developed  an  analytical  model  to  predict 
ths  fraction  of  tube  surface  that  is  flooded  luring  conden- 
sation on  a  horizontal,  integral-fin  tube.  They  found  that 
the  vertical-rise  height  of  condensate  on  a  finned  tube  was 
the  same  as  that  obtained  on  a  vertically  oriented,  flat 
finned  plate  that  was  obtained  by  splitting  and  unrolling  an 
identical  finned  tube.  Based  on  this  observation,  they  used 
capillary  equations  that  predict  the  liquid  rise  on  a 
vertical  U-shaped  channel,  and  assuming  negligible  vapor 
shear,    they    predicted    the    condensate    retention    angle    to    be: 
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where 


e   = 


t    = 


9   = 


s    = 


fin    heig  lit, 

fin    thickness, 

density    of    condensate, 

acceleration    of    gravity, 

fin    spacing,    and 

root    diameter  of   tube. 


This  result  shows  that  the  condensate  retention  angle 
increases  both  with  increasing  fin  density,  and  with 
increasing    surface    tension-to-density    ratio.  It    should    be 

kept       in    mind,  however,       that      the      model    was       based    on      a 

vertical  surface,  so  it  should  only  be  used  for  angles  less 
than  about  30  degrees.  For  angles  below  this,  experimental 
results  involving  the  use  of  water,  R-11,  P-12,  ammonia,  and 
n-pentane  were  predicted  to  within  10  percent.  It  is  worth 
mentioning  that  equation  (2.7)  had  been  derived  in  1982, 
according    to   Russian    literature,    by    Rifert    [  12]. 

Like  Rudy  and  Webb  [10],  Owen  et  al.  [13]  also  recog- 
nized the  need  to  consider  condensate  retention  while 
analyzing  finned  tubes.  In  order  to  correct  for  this,  an 
assumption  was  made  that  the  condensate  retention  angle  was 
independent  of  condensation  rate,  so  a  static  analysis  was 
performed.         A    simple   force      balance    between    surface    tension 


27 


an3  gravitational   forces  resulted   in  an   equation  for   the 
condensate  retention  angle  as  shown  below: 


ros 
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vhere 


Df  =  overall  diameter  of  fins. 


This  equation  is  the  same  as  equation  (2.7),  except  that 
equation  (2.8)  is  independent  of  fin  thickness  (t) . 

In  1983,  Honda  et  al.  [3]  performed  a  theoretical  anal- 
ysis to  determine  the  condensate  retention  angle.  Using  an 
iterative  numerical  scheme,  they  found  the  solution  to  agree 
with  equation  (2.8)  obtained  earlier  by  Owen  et  al. .  Using 
data  of  their  own,  Honda  et  al.  verified  a  close  agreement 
between  the  predicted  and  experimental  values  of  the  conden- 
sate retention  angle. 

Yau  et  al.  [  1 ]  conducted  experiments  to  determine  the 
effects  of  removing  retained  condensate  from  finned  tubes  by 
installing  thin  metal  drainage  strips  attached  edgewise  to 
the  bottom  tube  surface.  Experiments  were  conducted  where 
the  vapor  velocity  and  fin  pitch  were  varied  on  finned  tubes 
with  and  without  drainage  strips.  They  showed  a  significant 
reSuction  in  the  in  the  condensate  retention  angle  when  the 
finned  tubes  were  fitted  with  drainage  strips.  Condensate 
retention  angles  for  condensation  of  steam,  ethylene  glycol, 
ani  Fc-1-13  on  finned  tubes  with  drainage  strips  were  found  to 
fit  the  empirical  relation  listed  below  to  within  ±  10 
percent: 


28 


i'    =    cos 


-1 


1.66  o 


(2.9) 


In  1985,  Rudy  and  Webb  [14]  expanded  their  1983  model  in 
order  to  predict  the  condensate  retention  angle  for  hori- 
zontal finned  tubes  with  fins  of  arbitrary  shape.  As 
before,  their  model  was  based  on  capillary  equations  that 
predict  the  amount  of  liquid  rise,  and  negligible  vapor 
shear.  The  resulting  equation  for  the  condensate  retention 
angle  is  given  by: 
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where 


P^  =  wetted  perimeter  of  fin  cross  section, 

t.  =  fin  base  thickness,  and 

A   =  profile  area  of  fin  over  fin  cross  section 


In  order  to  test  the  validity  of  equation  (2.10),  it  was 
compared  to  experimental  data  for  horizontal  finned  tubes 
with  fins  of  primarily  trapezoidal  cross  section,  and  fin 
densities  ranging  from  630  to  1614  fins/m.  Steam,  F-11,  and 
n-pentane  were  used   as  the   condensing  fluids.     Equation 

(2.10)  was  shown  to  predict  the  experimental  data  to  within 
±  13  percent.   Note  that  equation  (2.10)  reduces  to  equation 

(2.8)  for  fins  of  rectangular  cross  section. 
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C.       THEORETICAL    MODELS 

In  1947,  Beatty  and  Katz  [15]  performed  a  number  of 
experiments  where  various  refrigerants  were  condensed  on 
single  finned  tubes.  To  predict  the  film  coefficients,  they 
started  with  the  Nusselt  equations  for  condensation  on  a 
horizontal  tube  and  on  a  vertical  plate,  and  with  the 
respective  proportions  of  horizontal  tube  area  (between 
fins)  and  vertical  fin  area,  they  calculated  the  total  heat 
transfer    rate.  From  this,    an   average    heat-transfer    coeffi- 

cient was  obtained  based  on  an  equivalent  tube  diameter. 
Their  final  expression  is  simply  the  Nusselt  equation  for  a 
smooth  horizontal  tube,  but  with  the  tube  diameter  replaced 
with  the  equivalent  tube  diameter.  The  leading  coefficient 
was    modified   to    fit    their   experimental    data    as    shown    below: 
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where 


AQ     =    surface    area    of    smooth    tube, 
h    =    average    vapor-side    heat-transfer   coefficient, 
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kf  =  thermal    conductivity    of   condensate, 

Pv  =  density    of    vapor, 

hf  =  specific    enthalpy    of    vaporization, 

yf  =  viscosity   of    condensate, 

At  =  vapor-side    temperature    drop, 

De  =  equivalent    tube   diameter, 

n  =  fin   efficiency, 

Ag  =  surface   area      of    smooth    tube, 

Af  =  total    surface  area   of    finned    tube,    and 

A ef f  =  effective  area   of    finned    tube. 

Accuracy  of  this  equation  was  claimed  to  be  better  than  1 1 
percent  for  a  wide  variety  of  nonagueous  fluids.  It  should 
be  noted,  however,  that  this  model  ignores  the  effect  of 
surface-tension  forces;  so,  it  is  valid  only  for  low- 
surface-tension    fluids   such    as   refrigerants. 

Analytical  and  experimental  studies  of  condensation  on 
horizontal  tubes  with  trapezoidally-shaped  fins  were 
performed  by  Zozulya,  Karkhu,  and  Borovkov  [16,17]  some 
years    later    in    the    1970s.  Their   experiments   confirmed    the 

need  to  consider  surface-tension  forces  in  addition  to 
gravitational  forces  when  developing  mathematical  models. 
The  analytical  solutions  were  based  on  the  following  assump- 
tions: 1)  the  thin  condensate  film  on  the  fins  was  treated 
as  a  laminar  boundary  layer  with  a  pressure  gradient  along 
the  fin  profile  caused  by  surface-tension  forces;  2)  the 
effect  of  gravitational  and  inertial  forces  on  the  motion  of 
the      film  along      the      side    surfaces      of      the      fins    into      the 
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condensate-filled  trough  was  neglected;  3)  the  motion  of 
condensate  in  the  tiough  area  is  laminar  and  produced  by 
gravity;  and  4)  no  condensation  takes  place  on  the  flooded 
portion  of  the  tube.  Differential  equations  were  obtained 
for  the  height  of  condensate  between  fin  bases,  and  tempera- 
ture distribution  along  the  fin  height.  With  the  aid  of  a 
computer,  solutions  for  these  equations  were  obtained  using 
numerical  methods,  leading  to  the  expressions  below: 
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where 


G  =  condensate  flow  rate, 

tt  =  fin  tip  width, 

a  =  fin  semivertex  angle,  and 

Z  =  dimensionless  depth  of  condensate 
between  fins. 

Although  the  above  eguations  appear  to  be  fairly  simple,   it 
must  be  kept  in  mind   that  rather  complex  computer  solutions 
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are  required  to  obtain  values  for  Z.  To  check  the  validi+y 
of  equation  (2.15),  (2.16)  and  (2.17;  ,  experiments  were 
performed  with  four  different  finned  tubes  to  condense  both 
steam  and  Ft-113  vapor.  Agreement  to  within  5  percent  was 
reported. 

In  1977,  Nader  [18]  presented  an  analytical  model  to 
predict  condensation  on  a  vertical,  single  fin  attached  to  a 
horizontal  tube.  This  model  differed  from  those  discussed 
above  in  that  the  fin  temperature  was  allowed  to  vary  along 
the  height  of  the  fin.  For  the  special  case  where  the  fin 
temperature  is  constant,  the  Nusselt  eguation  was  obtained. 
A  "condensation  efficiency"  was  then  proposed  to  account  for 
a  variable  fin  temperature,  and  was  defined  as: 
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where 


C,  =  constant  of  proportionality,  and 

k   =  thermal  conductivity  of  fin. 

m 


With  this  definition,  the  rate  of  condensate  formation  and 
the  rate  of  heat  transfer  (with  variable  fin  temperature) 
were  found  to  be: 
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Gf   =  rate  of  condensate  formation,  and 
Q  =  total  heat-transfer  rate. 


Thus,  with  a  table  of  values  for  F.  and  F  obtained  numeri- 
cally, the  rate  of  heat  transfer  and  condensate  formation 
could  be  calculated  by  a  fairly  simple  expression  provided 
the  condensation  efficiency  is  greater  than  0.8.  However, 
for  values  below  0.8,  the  accuracy  of  equations  (2.21)  and 
(2.22)  diminishes  quite  rapidly. 

In  1979,  Patankar  and  Sparrow  [19]  treated  the  case  of  a 
fin  attached  to  a  cooled  vertical  plate  or  a  cylinder  as  a 
three-dimensional  problem.  A  thin  fin  was  assumed,  and 
temperature  variations  across  the  thickness  were  ignored  as 
was  lengthwise  conduction.  Additionally,  all  heat  transfer 
was  assumed  due  to  condensation  only.  Compared  to  the  ideal 
fin  solution,  results  obtained  by  Patankar  and  Sparrow  indi- 
cated that  a  significant  error  existed  when  using  the  ideal 
fin  solution  on  all  but  very  long  fins.  Based  on  their 
results,  the  following  relations  were  proposed: 
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n  =  0.6171  hfr.  p>  pf  e 
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2    kf    el 


z  (2.  24) 


where 


g    =    heat    flux,    and 
z    =    axial   coordinate. 


According  to  Patankar  and  Sparrow,  these  eguations  should  be 
valid  for  most  practical  applications;  however,  no  experi- 
mental   data    were    available    at    the    time    to    confirm    this. 

In  1980,  Rifert  [20]  analyzed  condensation  of  vapor  on 
horizontal  finned  tubes  enhanced  by  the  effect  of  surface- 
tension  forces.  In  his  analysis,  he  divided  the  tube  into 
flooded  and  unflooded  zones,  and  solved  a  two-dimensional 
form    of   the      energy    eguation    for    each    zone.  The    mean    heat 

flux  was  then  determined  by  integrating  over  each  zone  and 
the  tube  length.  In  cases  where  condensate  is  retained  in 
more  than  half  of  the  tube  perimeter,  Rifert  points  out  that 
a  three-dimensional  form  of  the  energy  eguation  must  be 
used.  Solutions  to  these  eguations  revealed  that,  in  most 
cases,  the  fin  temperature  is  very  nonuniform,  so  the  mean 
integral  heat-transfer  coefficient  and  its  correlations 
should  not  be  used.  The  above-mentioned  zone-by-zone  anal- 
ysis   was,    therefore,    recommended. 

Based  on  their  study  of  condensate  retention  mentioned 
earlier    in    section    B,      Rudy    and      Webb   [10]    proposed    that    the 
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Beatty      and      Katz    [15]      model      be      modified   to      account      for 
condensate   retention,    as    shown    below: 


h    =  b 


PK 


ir    _    ty 


(2.25) 


where  li    is  computed  using  equation  (2.11).    This  equation 

BK 

neglects  any  heat  transfer  through  the  flooded  portion  of 
the  tube.  As  a  result,  it  was  shown  to  underpredict  the 
average  heat-transfer  coefficient  of  condensing  R-11  by  as 
much  as  30  percent  when  a  significant  amount  of  retained 
condensate  was    present. 

Using  the  condensate  retention  angle  they  had  developed 
earlier  (see  section  B)  ,  Owen  et  al.  [13]  divided  a  hori- 
zontal finned  tube  into  an  upper  unflooded  section,  and  a- 
lower  flooded  section.  Unlike  Rudy  and  Febb  [10],  however, 
they  accounted  for  heat-transfer  through  the  flooded  portion 
of  the  tube  by  noting  that  condensation  occurs  on  both  the 
surface  of  the  retained  condensate  and  the  fin  tips.  An 
effective  thermal  resistance  was  obtained  over  this  portion 
of  the  tube  by  assuming  parallel  paths  for  heat-transfer 
through   the      fins   and   retained   condensate.  The    respective 

thermal  resistances,  which  were  added  in  parallel,  were  then 
added  in  series  with  the  tube  wall  thermal  resistance.  The 
resulting   heat-transfer   coefficient    was    then    found   to   be: 


h   = 


(tt  -    i^ 


+   *h 


(2.26) 
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where    hy    is   computed    using    equation     (2.11)  #    and 


Leff 


+  i.^7Q    I  vr   ^  n" 


1/4-1 


-l 


(2.27) 


k    ...     =    effective   thermal   conductivity 


eff 


This    equation      was    found    to      precict    all    the      available    lata 
with    an    accuracy    of    better    than    30    percent. 

Fudy  and  Webb  [21]  used  equation  (2.7)  a  short  time 
later  while  developing  an  equation  for  the  heat-transfer 
coefficient.  The  model  for  this  undertaking  was  based  on 
surf ace-tension- driven  radial  flow  of  condensate  on  the  fin 
surface,  with  a  linear  pressure  gradient  along  the  fin 
height,  and  gravity-drained  flow  of  condensate  in  the 
channel  between  fins.  The  Nusselt  equation  for  horizontal 
tubes  was  used  for  the  tube  area  between  fins,  while  the  fin 
surface  was  treated  by  an  equivalent  gravity  model  developed 
by  Rebb  et  al.  [22]  and  R  u3  y  [23]  earlier.  As  before,  vapor 
shear  and  heat  transfer  through  the  flooded  portion  were 
neglected,    yielding    a  heat-transfer    coefficient    of: 
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where 


I  =  length  of  tube, 
Abt  =  surface  area  of  tube  between  fins,  and 
A f t  =  fin  surface  area. 

This  expression  provided  an  accuracy  of  better  than  1 0 55  for 
condensation  of  R-11  on  short,  finely-spaced  fins,  but  the 
accuracy  dropped  sharply  for  larger  fins  spaced  further 
apart.  This  was  apparently  due  to  the  assumed  linear  pres- 
sure gradient  on  the  fin  surface,  and  breakdown  of  the  model 
as  gravity  forces  became  dominant.  The  use  of  eguation 
(2.28)  should  be  restricted  to  fin  densities  from  1200  to 
1400  fins/m,  and  fin  heights  of  less  than  1  mm. 

In  1984,  Honda  et  al.  [24]  developed  a  model  to  predict 
the  heat-transfer  coefficient  that  took  surface-tension 
effects  and  the  n on- isothermal  behavior  of  fins  into 
account.  In  order  to  obtain  a  solution,  however,  iterative 
finite  difference  techniques  that  required  a  considerable 
amount  of  computer  time  were  used.  Nonetheless,  they  showed 
that  this  model  was  able  to  predict  the  average  heat- 
transfer  coefficient  for  most  available  experimental  data 
(which  included  11  fluids,  and  22  tubes)  within  20  percent. 
Predictions  involving  the  condensation  of  steam  provided  the 
largest  errors  (as  much  as  40/*)  ,  due  in  part  to  the  high 
surface  tension  of  water.  This  reflects  the  fact  that 
surface-tension  effects  are  still  not  completely  accounted 
for.  Compared  to  the  other  models  discussed  thus  far, 
however,  the  model  of  Honda  et  al.  shows  considerable 
promise  in  spite  of  its  complexity. 

In  1983,  Adamek  [27]  defined  a  family  of  condensate 
surface  profiles  whose  curvatures  are  given  by: 
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U  +   1) 


(2.29) 


where 


m 


m 


S  = 


rotation  angle  of  normal  to  condensate  film  surface, 

length  of  convex  surface  of  fin  condensate  film, 

coordinate  along  condensate  film  surface  from 
fin  tip. 


The  parameter  e,  in  eguation  (2.29)  characterizes  the  aspect 
ratio  of  the  fin  cross  section.  A  number  of  E  values  and 
their  corresponding  condensate  surface  profiles  are  shown  in 
Figure  2.3.  With  the  above  information,  Adamek  obtained  an 
expression  for  the  heat-transfer  coefficient  of  the  fin 
surface  as  shown  below: 


h   =  2.UH  ^ 
'  S 
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a  h 
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om  Sm  Pf  (€+  1) 


Uf  kf  AT 


a   +  2^ 


(2.30) 


These  fin  profiles  were  used  by  Webb  et  al.  [25]  and  Budy 
et  al.  [26]  in  1985  to  develop  a  model  for  condensation  on 
horizontal  finned  tubes.  This  moiel  divided  the  tube  into 
flooded,  and  unflooded  regions,  and  further  divided  the 
unflooded  region  into  tube  (area  between  fins)  and  fin 
areas.  To  predict  heat  transfer  from  the  unflooded  fin 
surfaces,  they  used  the  heat-transfer  coefficient  proposed 
by  Adamek.  The  Adamek  fin  profile  used  to  approximate  their 
trapezoidal   fins  is   shown   in   Figure  2.  4.     The   Nusselt 
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Figure  2.3    Adamek  [27]  Condensate  Surface  Profiles 
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Figure   2.4        Fin    Geometry    for   the    Webb   et   al.    Hodel  [22] 
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equation,  written  in  terms  of  the  film  Peynolds  number,  was 
used  to  predict  heat  transfer  from  the  unflooded  tube 
surface  between  fins.  Writing  the  Nusselt  equation  this  way 
takes  into  account  the  additional  condensate  thickness  that 
results  from  the  drainage  of  condensate  formed  on  the  fin 
surface.  For  the  flooded  region,  a  two-dimensional  computer 
code  was  used  to  solve  for  the  heat  flux  into  the  tube-side 
coolant  ( q  „) ,  which  was  compared  to  the  limiting  case  of 
heat    flux    (qbl)  assuming    zero    fin    thickness.  This    ratio 

(<t>  =  c3b2//^bl  '  was  use(^  along  with  a  linear  temperature 
profile  across  the  condensate  film  to  establish  a  heat- 
transfer  coefficient  for  this  region.  The  resulting  average 
heat-transfer    coefficient    is  given    below: 


h       -  h  n     - 
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+    n   h 
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(2.  31) 


=    1.514 
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h.    = 


(2.33) 


n     -  1  -  (1  -  rO         Aff 

Aht   +  Afr 


(2.  34) 


where    h f  is   calculated   using    equation    (2.30),    and 


A        =    smooth    tube    area    based   on    the    fin    diameter, 

sf 


Re-    =    film   Reynolds   number, 


41 


<\>    =    ratio    of    tube  side    heat    flux    with    fins    to    tube 
side   heat    flux   with    fins    of    zero    thickness. 

A    detailed    step-by-step    procedure,    recommended    for    obtaining 
a   solution    of   equation    (2.3  1),    is    provided    in   Appendix    A. 

To  allow  for  heat  transfer  from  the  fin  tip,  and  to 
account  for  decreased  heat  transfer  due  to  the  thick  conden- 
sate film  in  the  trough  area,  Webb  et  al.  replaced  S^  with 
sm  +  t/2  -  6  ,  where  6r  is  the  average  condensate  film  thick- 
ness   in    the   trough    area    in    the   unflooded   region: 


Since  the  effects  due  to  t/2  and  6r  in  the  modified  Sm  tend 
to  cancel  each  other,  Webb  et  al.  reported  only  a  negligible 
change   in    Sm.  This   model    was      shown    to    predict      the    heat- 

transfer  coefficient  for  E-11  condensing  on  horizontal 
finned  tubes  with  fin  densities  of  748,  1024,  and  1378 
fins/m    within    20    percent. 
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III.     DESCRIPTION    OF    TEST    APPARATUS 

A.       TEST    APPARATUS 

The    same   test    apparatus    used      by    Georgiadis    [7]    was   used 
for    this   investigation.  A    schematic    of    this      apparatus   is 

shown  in  Figure  3.1.  Steam  was  generated  in  the  boiler, 
which  consisted  of  a  304.8  mm  (12  in.)  Pyrex  glass  section 
containing  ten  4000-Watt,  480-Volt  Watlow  immersion  heaters. 
After      passing      through      a    304-8    mm       (12    in.)  to       152.4    mm 

(6  in.)  reducing  section,  the  steam  flowed  upward  through  a 
2.44  m  (8  ft.)  long  section  of  Pyrex  glass  piping.  At  this 
point,  a  180-degree  bend  in  the  piping  re-directed  the  steam 
downward,  where  after  1.52  m  (5  ft.)  of  Pyrex  glass  piping, 
it  entered  the  stainless  steel  test  section  illustrated  by 
Figure  3.2.  The  test  tube  was  mounted  horizontally  in  the 
test  section  behind  a  view  port,  which  allowed  visual  obser- 
vation of  the  condensation  process.  Steam  not  condensed  by 
the  test  tube  continued  downward  to  the  auxiliary  condenser, 
where    it      was   condensed      by    two      9.5    mm     (3/8   in.)  diameter 

water-cooled  copper  tubes  helically  coiled  to  a  height  of 
457  mm  (18  in.).  The  condensate  was  then  returned  to  the 
boiler   by   gravity. 

Filtered    tap    water      was    used   to    cool    the   test      tube    on   a 
once-through   basis.  This    water   was      first    collected      in    a 

large  sump  with  a  capacity  of  about  0.4  cubic  meters  (Figure 
3.3),  then  pumped  through  a  flow  meter,  and  the  test  tube  by 
two  centrifugal  pumps  connected  in  series.  A  valve  on  the 
discharge  side  of  the  second  pump  allowed  the  velocity  of 
water  flowing  through  the  test  tube  to  be  varied  from  0  to 
4.4  m/s  (14.4  ft/sec).  The  auxiliary  condenser  was  cooled 
by  a  continuous  supply  of  tap  water,  which  was  throttled  to 
control   the    internal    pressure    of    the    test    apparatus. 
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The  system  used  to  remove  non-condensing  gases  is  shown 
in  Figure  3.3.  The  vacuum  pump  continually  drew  a  sample 
mixture  of  any  possible  non-condensing  gases  and  moisture 
from  the  auxiliary  condenser.  This  mixture  was  then  passed 
through  a  heat  exchanger,  where  the  moisture  was  condensed 
before  it  entered  the  sump,  using  the  same  filtered  tap 
water  described  earlier. 

B.  INSTBUMENTATION 

The  electrical  power  input  to  the  boiler  immersion 
heaters  was  controlled  by  a  panel-mounted  potentiometer.  A 
converter  with  an  input  voltage  of  440  VAC  generated  a 
signal  which  was  fed  to  the  data  acquisition  system  to 
calculate  the  boiler  input  power.  A  more  detailed  descrip- 
tion   of   the  boiler   power   supply  is   provided   by   Poole   [6]. 

The  internal  system  pressure  was  measured  manually  using 
a  0-tube,  mercury-in-glass  manometer  graduated  in  millime- 
ters. Steam,  condensate  return,  and  ambient  temperatures 
were  measured  using  calibrated  copper-constantan  thermocou- 
ples. These  thermocouples  had  an  accuracy  within  ±  0.1  K. 
The  cooling  water  temperature  rise  was  measured  by  two 
Hewlett-Packard  (HP)  2804A  quartz  thermometers,  along  with  a 
10- junction,  series-connected  copper-constantan  thermopile 
as  a  backup.  Throughout  this  investigation,  the  quartz 
thermometers   and   the   thermopile   agreed   to   within    ±    0.03   K. 

C.  DATA  ACQUISITION  SYSTEM 

To  monitor  the  system  temperatures  and  boiler  input 
power  (using  the  converter  signal),  an  HP  9826A  computer  was 
used  to  control  an  HP  3497A  Data  Acquisition  System.  Raw 
data  were  processed  immediately  and,  at  the  same  time, 
stored  on  diskett  for  reprocessing  at  a  later  time. 
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D.  TUBES   TESTED 

For  this  thesis  effort,  eighteen  copper  tubes,  and  a 
stainless  steel  tube  were  manufactured.  They  are  presented 
in  Table  I.  The  first  eight  tubes  with  rectangularly-shaped 
fins  were  tested  in  order  to  complete  the  sequence  of  data 
initiated  by  Seorgiadis  [7].  This  completed  the  systematic 
variation  of  fin  spacing,  thickness,  and  height  so  that  the 
optimum  dimensions  could  be  determined.  Figure  3.4  shows  a 
photograph  of  the  four  tubes  with  fin  thickness  of  1.0  mm, 
fin  height  of  0.5  mm,  and  spacings  of  1.0,  1.5,  2.0,  and  4.0 
mm    (tube   numbers   26    to  29).  Also    shown    is   the    insert   used 

to  enhance  the  inside  heat- transfer  coefficient,  which  will 
be  discussed  in  more  detail  in  the  next  chapter.  Figure 
3.5(a)  shows  a  cross  section  of  a  rectangular  fin  with  a  fin 
height  of  1.5  mm,  a  fin  spacing  of  1.5  mm,  and  a  fin  thick- 
ness of  1.0  mm.  In  addition  to  the  tubes  with 
rectangularly-shaped  fins,  six  tubes  were  tested  to  study 
the  effect  of  fin  shape  on  the  heat-transfer  performance. 
Tubes  number  40  and  42  were  made  of  stainless  steel  to  test 
the  effect  of  fin-metal  thermal  conductivity  on  the  heat- 
transfer  performance.  Due  to  limited  machine-shop  facili- 
ties available,  the  "parabolic"  fins  on  two  of  the  tubes 
were  not  truly  parabolic,  as  shown  in  Figure  3.5(b).  A 
cross  section  of  a  commercially-available  finned  tube 
produced  by  the  Wolverine  Division  of  United  Oil  Products  is 
also  shown  in  Figure  3.5(c)  for  comparison.  The  smooth 
tubes  were  tested  in  order  to  determine  the  water-side  heat- 
transfer  coefficients,  and  to  serve  as  a  comparison  for  the 
finned    tubes. 

E.  VACUUM    INTEGRITY 

For  any  condensation  experiment,  especially  when  oper- 
ating  under    vacuum    conditions,      a      leak-free   apparatus   is  of 
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Figure  3.5    Cross- Sectional  Photographs  of  (a).  Rectangular 
Fins,  (b)  "Parabolic"  Fins,  and  (c)  "Wolverine"  Fins. 
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vital  importance.  The  presence  of  any  leaks  would  allow  the 
admission  of  non-condensing  gases,  whose  deleterious  effects 
would  invalidate  any  data  being  collected.  To  prevent  this 
problem,  vacuum  integrity  tests  were  performed  on  the  test 
apparatus  both  before  and  after  the  data  acquisition  phase 
of  this  thesis  effort  (a  time  period  of  about  four  months) . 
The  first  test  revealed  a  leak  rate  represented  by  a  pres- 
sure rise  of  1.0  mmHg  over  a  period  of  five  days,  while  the 
pressure  rise  during  the  final  test  was  1.0  mmHg  over  a 
24-hour  period.  These  tests  were  performed  at  an  absolute 
pressure      of   about      85  mmHg.  This      negligible    leak      rate, 

together  with  continuous  venting  of  the  test  apparatus 
during  data  runs,  eliminated  the  deleterious  effects  of 
non-condensing  gases  on  the  heat-transfer  results  reported 
in    this    thesis. 
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17.  DATA  COLLECTION  AND  REDOCTION 

A.   SYSTEM  OPERATION 

In  order  to  ensure  filmwise  condensation,  all  tubes  were 
treated  with  a  solution  containing  equal  parts  of  sodium 
hydroxide  and  ethyl  alcohol.  The  tubes  were  placed  in  a 
steam  bath  and  several  layers  of  the  solution  were  applied 
with  a  tooth  brush  every  10  minutes.  After  about  one  hour, 
the  tube  surfaces  were  completely  blackened.  The  same 
procedure  was  repeated  prior  to  each  data  run,  but  the  tube 
was  only  steamed  for  ten  to  fifteen  minutes.  Following  the 
steam  bath,  the  tubes  were  rinsed  with  tap  water  before 
installing  them  in  the  test  apparatus.  Georgiadis  [7] 
discusses  this  procedure  in  more  detail. 

The  test  apparatus  was  brought  to  operating  pressure  and 
temperature  using  the  procedures  provided  by  Georgiadis  [7], 
Steady-state  conditions  were  assumed  to  occur  once  the  test 
section  steam  temperatures,  as  well  as  the  cooling  water 
inlet  and  outlet  temperatures  stabilized.  Two  sets  of  data 
were  then  taken  for  each  of  the  heat  fluxes,  which  were 
determined  by  the  flow  rate  of  cooling  water  through  the 
test  tube.  Starting  at  a  flow  rate  of  80  percent  (4.4a  m/s 
for  19  mm  O.D.  tubes,  and  3.92  m/s  for  13.7  and  15.8  mm  C.D. 
tubes),  flow  rates  were  reduced  to  70,  60,  45,  35,  26,  an  1 
20  percent,  then  brought  back  up  to  55  and  90  percent. 
These  settings  were  selected  since  they  provided  nearly 
equally-spaced  heat  flux  values.  It  should  be  noted  that 
continuous  adjustments  of  cooling  water  flow  through  the 
auxiliary  condenser  were  required  to  maintain  system  pres- 
sure as  flow  rates  through  the  test  tube  were  changed. 
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Visual  observations  of  condensation  taking  place  on  the 
test  tube  were  made  on  a  regular  basis  to  ensure  uniform 
filmwise  condensation.  A  more  reliable  indication  of  this, 
however,  was  the  trend  of  data  taken  as  the  pooling  water 
flow    rates      were   increased   back  up      to   80    percent.  If  the 

final  value  of  the  steam-side  heat- transfer  coefficient  at 
80  percent  was  significantly  higher  than  those  at  the  begin- 
ning of  the  data  run,  dropwise  condensation  was  assumed  to 
have    occurred,         and    the   data      were    disregarded.  All   data 

presented  in  this  thesis  displayed  less  than  5%  disagreement 
in  the  steam-side  heat-transfer  coefficient  between  initial 
and    final    data    sets. 

B.  THE    DROPWISE    CONDENSATION    PROBLEM 

The  purpose  of  this  thesis  was  to  take  strictly  filmwise 
condensation  3ata,  using  primarily  copper  tubes.  Due  to  the 
poor  wettability  of  copper  with  water,  especially  when  even 
minute  amounts  of  contamination  are  present,  there  is  a 
tendency  to  condense  steam  in  the  partial  dropwise  mode. 
Dropwise  condensation  is  far  more  effective  than  filmwise 
condensation;  so,  its  presence  can  lead  to  large  errors  in 
the  data. 

By  following  procedures  set  forth  by  Georgiadis  [7],  the 
dropwise  problem  was  minimized.  During  the  initial  testing 
phase  of  this  thesis,  however,  a  small  amount  of  dropwise 
condensation  was  encountered.  This  problem  was  soon  traced 
to  the  boiler  water,  which  had  not  been  changed  for  some 
time.  By  changing  the  water  after  every  three  or  four  runs, 
the  problem  was  eliminated. 

C.  STEAM  VELOCITY  LIMITATIONS 

A  major  assumption  in  the  Nusselt  theory  is  that  vapor 
velocity   does      not    induce   shear      forces   on      the    vapor-licuid 
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interface.  Since  system  pressure  in  the  test  apparatus  was 
controlled  by  condensing  steam  in  the  auxiliary  condenser,  a 
reasonable  steam  velocity  was  always  present.  Since  steair 
velocity  (i.e.,  steam  flow  rate)  depends  on  boiler  power 
input,  decreasing  the  power  results  in  lower  vapor  veloci- 
ties. Under  vacuum  conditions,  (about  85  mmHg)  a  maxiDum 
velocity  of  about  8  m/s  was  possible  at  full  boiler  power 
(36  kW) .  However,  at  these  high  velocities,  the  entire 
apparatus  experienced  considerable  vibration.  As  discussed 
by  Georgiadis  [7],  a  steam  velocity  of  2  m/s  under  vacuum 
provided  the  most  stable  operating  conditions.  This  corre- 
sponding boiler  power  (about  9  kW)  was  low  enough  to  mini- 
mize system  vibrations,  yet  high  enough  to  avoid  rapid 
fluctuations  in  condenser  pressure  due  to  the  intermittent 
break-up  of  vapor  bubbles  in  the  boiler.  To  allow  a  suffi- 
cient amount  of  steam  to  reach  the  auxiliary  condenser  when 
operating  at  atmospheric  pressure,  the  maximum  boiler  power 
was  required.  This  resulted  in  a  steam  velocity  of  approxi- 
mately 1  m/s. 

D.   DATA  REDUCTION 

All  the  programs,  property  functions,  and  calibration 
data  used  for  this  thesis  effort  were  essentialy  the  same  as 
those  used  by  Georgiadis  [7]  and  Poole  [6].  Because  six  of 
the  tubes  (tube  numbers  39  thru  44)  tested  during  this 
effort  had  inside  diameters  different  than  any  tube  tested 
previously,  and  two  (tubes  number  40,  and  42)  were  manufac- 
tured from  stainless  steel,  some  minor  modifications  were 
necessary  for  the  data  reduction  program  (program  DRP6) . 
These  modifications  included  options  for  different  tube 
diameters,  tube  inserts,  thermal  conductivity,  mixing- 
chamber  calibrations,  and  leading  constants  for  the  inside 
heat-transfer  coefficients. 
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The  mixing- chamber  (see  Figure  3.2)  calibration  was 
required  to  account  for  the  temperature  rise  that  occurred 
across  the  mixing  chamber,  due  primarily  to  turbulent 
viscous  dissipation.  The  purpose  of  the  mixing  chamber  is 
to  ensure  a  uniform  temperature  at  the  crooling  water 
discharge,  so  that  accurate  temperature  measurements  can  be 
male.  Since  tubes  with  different  inside  diameters, 
requiring  different  inserts,  result  in  different  pressure 
drops,  calibrations  had  to  be  performed  for  each  inside 
diameter. 

A  modifiel  Wilson  Plot  program  [7]  was  used  to  compute 
the  leading  constants  for  the  Sieder-Tate  equation  used  to 
determine  the  water-side  heat-transfer  coefficient.  Data 
taken  on  smooth  tubes  were  used  to  determine  the  leading 
constants,  which  were  then  used  to  reduce  data  taken  on 
similar  tubes  with  identical  internal  dimensions. 
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V.  REVOLTS  AND  DISCISSION 

A.  INTRODUCTION 

During  this  thesis  effort,  a  number  of  data  runs  were 
made  using  the  procedures  describe!  in  Chapter  IV.  For  both 
vacuum  and  atmospheric  conditions,  every  tube  was  tested  at 
least  three  times,  on  different  days,  to  ensure  the  data 
were  repeatable.  Complete  filmwise  condensation  was  main- 
tained, and  the  computed  non-condensing  gas  concentration 
was  held  below  1.0  percent  for  all  accepted  data  runs.  It 
was  estimated  that  an  error  of  ±  0.2  K  in  the  steam  tempera- 
ture, or  an  error  of  ±  1.0  mmKg  in  the  system  pressure  would 
lead  to  an  error  of  ±  1.2  percent  in  the  non-condensing  gas 
concentration.  This  shows  that  the  computed  non-condensing 
gas  concentration  was  zero  to  within  the  accuracy  of  temper- 
ature and  pressure  measurements.  As  discussed  earlier  in 
Chapter  III,  the  test  apparatus  would  allow  only  a  negli- 
gible amount  of  non-condensing  gas  to  be  present  in  the 
apparatus.  The  non-condensing  gas  concentration  was 
computed  for  the  system  temperature  and  pressure  during 
every  data  run  to  ensure  no  major  leaks  developed.  A 
summary  of  tubes  tested  between  Georgiadis  [7]  and  this 
thesis  effort  and  the  resulting  enhancements  are  provided  in 
Table  II. 

B.  flATEB-SIDE    HEAT-TEANSFER    COEFFICIENTS 

The  water-side  heat-transfer  coefficients,  enhanced  by 
the  use  of  an  insert,  were  obtained  by  Georgiadis  [7]  for 
the  tubes  with  an  inside  diameter  of  12.7  mm.  Both  a 
"direct  method,"  and  a  "modified  Wilson  method"  were  used, 
and    are    discussed   below. 
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1 •       The    Direct    Method 

For  the  "direct  method,"  a  smooth  tube  was  manufac- 
tured vith  six  wall  thermocouples  inserted  around  its  perim- 
eter. During  data  runs  under  vacuum  and  atmospheric 
conditions,  readings  from  these  six  thermocouples  were 
obtained  and  averaged,  in  order  to  determine  an  average  tube 
wall  temperature  for  each  heat  flux.  It  should  be  noted 
that  the  wall  temperature  was  quite  non-uniform  around  the 
tube  perimeter.  Georgiadis  [7]  reported  up  to  an  18-K  drop 
in  the  wall  temperature  from  the  top  to  the  bottom  of  the 
tube.  Evaluation  of  the  uncertainty  resulting  from  this 
temperature  variation  was  considered  beyond  the  scope  of  the 
investigation.  The  average  wall  temperature  was  then  used 
along  with  the  cooling  water  inlet  and  outlet  temperatures 
to  determine  the  log-mean-temperature  difference  (LMTD) . 
With  values  for  heat  flux  and  LMTD  known,  the  value  of  the 
water-side  heat-transfer  coefficient  was  obtained.  This,  in 
turn,  was  used  to  determine  the  leading  coefficient  for  the 
Sieder-Tate-type    equation   shown   below: 


Nil    =    C   Re0'8    Prl/3 


0.  14 


+    B 


(5.1) 


A  more  detailed  step-wise  solution  procedure  is  given  by 
Georgiadis      [7].  For   the      previously-mentioned      tubes     of 

15.8  mm  and  19.05  mm  0. D. ,  Georgiadis  determined  the  coeffi- 
cient C  to  be  0.0635.  The  constant  B  =  26. H  is  an  addi- 
tional parameter  found  to  improve  the  fit  of  equation  (5. 1) . 
ITote  that  C  =  0.0635  is  about  2.5  times  greater  than  the 
well-known  Sieder-Tate  coefficient  of  0.027  [28]  for  long 
tubes  with  smooth  inner  walls.  This  is  due  primarily  to  the 
insert  used  to  enhance  the  water-side  heat-transfer 
coefficient. 
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2.       The    Modified    Wilson    Method 

A  complete  description  of  this  method  is  provided  in 
Appendix  B.  Briefly,  this  method  uses  a  Sieder-?ate-t ype 
equation  for  the  water-side  heat-transfer  coefficient,  and 
an  equation  (discussed  below)  for  the  steam-side  heat- 
transfer  coefficient.  Both  equations  contain  leading  coef- 
ficients that  must  be  determined  iteratively  while  being 
used  in  a  linear  equation  to  generate  the  Wilson  plot.  The 
iterations  were  performed  on  the  HP  9826A  computer  mentioned 
previously   in   Chapter  III,    using    program    "WILS0N6." 

The  "modified  Wilson  method"  used  during  this  inves- 
tigation consisted  of  one  minor  modification  from  that  of 
Georgiadis.  For  this  thesis,  a  Fujii-type  [29]  equation  was 
used  instead  of  the  Nusselt-type  eguation  for  the  steam-side 
coefficient.  This  modification  was  made  to  account  for  the 
small  steam  velocity  (about  1  m/s)  that  was  present  during 
the  runs.  For  this  reason,  use  of  the  Fujii-type  equation 
appears  to  be  more  accurate  than  the  use  of  the  Nusselt-type 
equation.  This  modification  resulted  in  a  slightly  higher 
value  (up  to  3  percent)  for  the  leading  coefficient  C.  The 
program  "T-7ILSON6"  described  earlier,  allows  an  option  for 
selecting  either  the  Fujii-type  or  Nusselt-type  equation  for 
the   steam-side   coefficient. 

Using  a  Nusselt-type  equation  for  the  steam-side 
coefficient,  Georgiadis  [7]  found  the  leading  coefficient 
(for  equation  (5.1))  C  to  be  0.071,  with  the  B  value  set 
equal  to  zero.  This  C  value  is  about  10  percent  higher  than 
the  value  (C  =  0.0635)  obtained  by  the  "direct  method,"  but 
since  the  "direct  method"  is  generally  felt  to  be  more  reli- 
able, the  values  of  0.0635  and  26.4  were  used  for  the 
constants  C  and  B  respectively  for  this  thesis  effort  with 
tubes   of    12.7   mm  inside   diameter    (see   Table    I). 
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3-      Water-Side    Coefficients    For    Thin- Wall   Tubes 

Due  to  the  thin  (i.e.,  0.5  mm  thickness)  tube  walls 
for  tubes  39  thru  42,  it  was  not  possible  to  manufacture  an 
instrumented  tube.  For  this  reason,  the  "modified  Wilson 
method,"  together  with  the  Fujii-type  equation  for  the 
steam-side  coefficient,  was  used  to  determine  the  leading 
coefficient  for  equation  (5.1)  for  both  the  copper  tube 
(tube  number  41)  ,  and  the  stainless  steel  tube  (tube  number 
42).  For   the      copper      tube,      a      value      of   C   =    0.0756      was 

obtained,  while  for  the  stainless  steel  tube,  a  value  of 
C  =  0.0688  was  obtained.  Although  the  inner  tube  diameters 
for  both  tubes  are  the  same  (i.e.,  13.5  mm),  the  much  higher 
thermal  conductivity  of  copper,  which  results  in  circumfer- 
ential heat  conduction  through  the  tube  wall,  leads  to  its 
higher  value  for  C.  A  summary  of  the  leading  coefficients 
used  in  equation  (5.1)  for  all  the  tubes  is  presented  in 
Table   III. 

C.       FEPEATABILITY   OF    DATA 

Since  the  reliability  of  data  taken  was  of  vital  impor- 
tance, all  data  runs  were  repeated,  as  Georgiadis  [7]  did, 
at  least  three  times  on  different  days.  The  computed  steam- 
side  coefficients  for  similar  conditions  (i.e.,  same  tube 
and  the  same  operating  conditions)  from  different  days 
agreed  to  within  ±  5  percent.  Additionally,  data  runs  were 
performed  on  two  finned  tubes  (tubes  number  6  and  17) ,  under 
similar  conditions  to  verify  the  repeatability  with  data 
taken  by  Georgiadis.  Figure  5.1  shows  that  the  experimental 
steam-side  heat-transfer  coefficients  of  Georgiadis  and 
those  obtained  during  this  investigation  agree  to  within  ±  5 
percent.  The  curves  shown  in  this  figure  (and  all  other 
figures  of  a  similar  nature  that  follow)  are  the  least- 
squares-fit   curves    according    to  the    following   equation: 
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q  =  a  ATb  (5>2) 


where  a  and  b  are  experimentally  determine!  constants. 
These  curves  were  selected  instead  of  those  based  on  the 
steam-side  coefficient  versus  heat  flux  since  the  uncer- 
tainty band  on  AT  is  more  uniform  than  the  uncertainty  band 
based  on  the  steam-side  coefficient.  A  brief  discussion  of 
the  uncertainty  bands  is  presented  in  Appendix  D,  along  with 
a  listing  of  the  uncertainty  analysis  program  "0NA6"  and  a 
few  sample  runs. 

D.   EFFECTS  OF  FIN  SPACING  AND  FIN  HEIGHT  ON  PEBFOBMANCE 

Data  were  taken  on  eight  tubes  (see  Table  II)  with  fins 
of  rectangular  shape  (tubes  number  26  thru  33)  under  vacuum 
(approximately  85  mmHg)  and  atmospheric  conditions.  Four 
tubes  had  a  fin  height  of  0.5  mm,  and  fin  spacings  of 
1-0,  1.5,  2.0,  and  4.0  mm,  while  the  remaining  four  tubes 
Lad  a  fin  height  of  1.5  mm,  and  the  same  sequence  of  fin 
spacings.  Georgiadis  [7]  reported  data  on  rec-tangular ly- 
shaped  fins  with  the  same  fin  spacings,  but  with  fin  heights 
of  1.0  mm,  and  2.0  mm.  A  complete  listing  of  tubes  tested 
by  Georgiadis  is  also  presented  in  Table  II  (tubes  1  thru 
25)  . 

Data  reductions  were  performed  on  the  KP  9926A  computer, 
using  a  program  named  "DRP6."  This  program  was  essentially 
the  same  one  used  by  Georgiadis  ("DRP4") ,  with  modifications 
to  allow  for  variations  in  tube  diameters,  and  tube 
material.  Data  obtained  for  all  tubes  tested  during  this 
thesis  effort  (tubes  number  26  thru  UU)  are  presented  in 
Appendix  C. 
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1  •      Effects    of    Fin   Spacing 

This  section  presents  results  showing  the  variations 
of  the  steam-side  heat-transfer  coefficient  with  heat  flux 
having    fin      spacing    as      a    parameter.  Figures    5.2      and    5.3 

present  data  for  tubes  with  a  fin  height  of  0.5  mm  under 
vacuum  and  at  atmospheric  pressure  respectively.  The  smooth 
tube  data  and  the  theoretical  prediction  of  Nusselt  are  are 
included   for   comparison.         In    a  similar   way,    Figures    5. 4   and 

5.5  present  data  for  tubes  with  a  fin  height  of  1.5  mm. 
These  figures  show  that  the  best  heat-transfer  performance 
occurs  for  the  tubes  with  a  fin  spacing  of  2.0  mm,  while  the 
worst  performance  occurs  for  the  tubes  with  a  fin  spacing  of 
4.0  mm.  Note  that  these  trends  are  the  same  both  under 
vacuum   and    at   atmospheric   conditions. 

The  comparison  of  finned  tubes  is  made  through  the 
enhancement  ratio.  This  ratio  is  defined  as  the  steam-side 
heat-transfer  coefficient  of  a  finned  tube  to  that  of  the 
smooth  tube  (same  diameter  as  the  finned  tube  root  diameter) 
at  the  same  heat  flux.  As  was  shown  by  Georgiadis  [7],  the 
uncertainty  in  calculating  the  steam-side  heat- transfer 
coefficient  increases  with  decreasing  heat  flux;  so,  heat 
flux  values  of  0.35  and  1.0  MW/m2  were  chosen  for  vacuum  "and 
atmospheric  conditions  respectively.  Cross  plots  of  the 
enhancement  ratio      versus    fin   spacing      are   shown      in   Figures 

5.6  and  5.7  for  fin  heights  of  0.5,  1.0,  1.5.  and  2.0  mm 
(curves  representing  e  =  1.0  and  2.0  mm  are  from  Georgiadis 
[7]).  Enhancements  as  high  as  4.8  and  6.4  can  be  seen  for 
vacuum  and  atmospheric  conditions  respectively.  Once  again, 
fin  spacings  of  approximately  2.0  mm  for  the  tubes  with  fin 
heights  of  0.5  and  1.5  mm  are  shown  to  provide  the  optimum 
performance.  Data  provided  by  Georgiadis  [7]  for  fin 
heights  of  1.0  and  2.0  mm,  however,  demonstrate  an  optimum 
heat-transfer   performance    for   fin    spacings   of    1.5   mm.         As   a 
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result,  there  does  not  appear  to  be  a  clear  trend  for  the 
optimum  fin  spacing  as  a  function  of  fin  height. 

On  one  hand,  the  difference  in  performance  between 
the  s  =  1.5  and  2.0  mm  finned  tubes,  as  well  as  the  s  =  0.5 
and  1.0  mm  finned  tubes  is  less  than  the  experimental  uncer- 
tainty (±  15  percent),  so  that  the  optimum  spacing  may  be 
somewhere  between  1.5  and  2.0  mm.  On  the  other  hand,  each 
data  point  in  these  figures  has  been  repeated  within  5 
percent  (i.e.,  computed  steam-side  heat-transfer  coeffi- 
cient) on  different  days.  Also,  as  discussed  earlier,  the 
data  taken  on  tubes  6  and  17  by  Georgiadis  were  repeated 
within  3  percent  during  this  investigation.  Therefore,  the 
slightly  disagreeing  trends  shown  by  Georgiadis  and  by  this 
investigator  (see  Figures  5.6  and  5.7)  may  be  justified  by 
the  very  complex  nature  of  the  condensation  phenomenon  on 
finned  tubes. 

The  presence  of  an  optimum  fin  spacing,  however,  is 
easy  to  understand.  The  1.0-mm  fin  spacing,  which  provides 
the  largest  increase  in  surface  area,  performed  worse  than 
the  2.0-mm  fin  spacing  due  to  the  relatively  large  conden- 
sate retention  angle  with  the  1.0  mm  fin  spacing.  Visual 
observations  showed  that  the  condensate  retention  angle  was 
around  110  degrees  for  the  tube  with  s  =  1.0  mm,  while  it 
was  about  40  degrees  for  the  tube  with  s  =  2.0  mm.  The 
additional  thermal  resistance  induced  by  this  layer  of 
condensate  overpowers  the  benefit  gained  from  increased 
surface  area,  so  the  heat-transfer  performance  is  reduced. 
As  fin  spacing  increases,  the  condensate  retention  angle 
decreases,  as  does  the  thickness  of  the  condensate  film 
between  fins  on  the  upper  portion  of  the  tube.  Beyond  fin 
spacings  of  2.0  mm,  the  decrease  in  surface  area  tends  to 
overshadow  the  other  effects,  resulting  in  decreased  heat- 
transfer  performance  compared  to  the  tube  with  the  "optimum" 
spacing. 
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2 .   Effects  of  Fin  Height 

As  mentioned  earlier,  fin  heights  of  0.5  and  1.5  mm 
were  examined  during  this  thesis  effort,  and  data  for  the 
tubes  with  fin  heights  of  1.0  and  2.0  mm  were  available  from 
Georgiadis  [7],  Table  II  shows  that  the  general  trend  is  an 
increase  of  enhancement  ratio  as  fin  height  increases,  which 
results  from  the  increase  of  surface  area  provided  by  the 
higher  fins. 

To  study  the  enhancement  obtained  beyond  the 
increase  in  surface  area  due  to  finning,  two  additional 
columns  are  provided  in  Table  II.  These  columns  represent 
enhancement  ratio/area  ratio  (Eo/Ar)  at  both  vacuum  and 
atmospheric  conditions.  As  can  be  seen,  Eo/Ar  values  range 
from  0.86  to  2.23  under  vacuum  conditions,  and  from  1.43  to 
3.20  at  atmospheric  pressure.  Figures  5.8  and  5.9  present 
plots  of  of  Eo/Ar  as  a  function  of  fin  spacing  under  vacuum 
conditions,  and  at  atmospheric  pressure  respectively.  The 
relatively  small  values  of  Eo/Ar  =  0.86  (vacuum)  and 
Eo/Ar  =  1.43  (atmospheric)  for  tube  number  15  are  due  to  the 
small  fin  spacing  which  results  in  a  completely  flooded 
tube.  On  the  other  hand,  the  values  of  Eo/Ar  =  2.23 
(vacuum)  for  tube  number  14,  and  Eo/Ar  =  3.20  (atmospheric) 
for  tube  number  28,  show  that  significant  enhancements 
beyond  the  area  ratio  are  obtainable  when  the  amount  of 
condensate  flooding  is  kept  small.  It  appears  that  surface- 
tension  forces  are  mostly  responsible  for  this  enhancement 
beyond  the  area  ratio. 

As  discussed  in  Chapter  II,  the  condensate  thinning 
provided  by  surface-tension  effects  is  related  to  the  pres- 
sure gradient  from  the  fin  tip  to  the  fin  root.  For  short 
fins,  this  pressure  gradient  is  relatively  larger  than  for 
high  fins,  resulting  in  improved  thinning  of  the  condensate 
film.    On  the  other  hand,  as  fin  height  increases,   surface 
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area  increases  and  the  condensate  retention  angle  decreases. 
Both  of  these  effects  lead  to  increased  heat-transfer 
perf ooance,  but  when  combined  with  the  poorer  thinning 
effect,  the  Eo/Ar  ratio  decreases  somewhat  for  the  larger 
fin    heights. 

E.       EFFECT    OF    FIN    GECHETRY    ON    PERFORMANCE 

In  order  to  study  the  effect  of  fin  geometry,  data  were 
taken  on  tubes  with  fins  of  triangular,  trapezoidal,  and 
parabolic  shapes  (tubes  number  34,  35,  37,  38)  ,  as  well  as 
on  a  spirally-finned  tube  with  threads  of  triangular  shape 
(tube      number    36).  Fin      dimensions       (i.e.,      fin      spacing, 

height,  and  thickness)  for  these  tubes  were  chosen  so  that 
the  resulting  total  surface  areas  would  be  approximately 
equal  to  that  of  the  "optimum"  tube  with  rectangularly- 
shaped  fins  and  a  fin  height  of  1.0  mm  (tube  number  6).  In 
this  manner,  the  effects  of  increased  surface  area  were 
eliminated,  allowing  a  comparison  between  tubes  where  fin 
geometry   was      the   only  variable.  A   commercially-available 

"Wolverine"  tube  (tube  number  43),  whose  fin  geometry  is 
shown  in  Figure  3.5(c),  was  also  tested.  Dimensions  for  the 
above   fins   are    given   in   Table   II. 

1 •      Iffect    of   Fin  Shape 

This    section    examines    the   effect      that    fin    shape    has 
on      the    heat-transfer      performance      of      finned    tubes.  The 

performance  of  tubes  number  34,  35,  36,  37,  and  43  under 
vacuum  conditions  is  shown  in  Figure  5.10,  while  Figure  5.11 
depicts      their    performance      at      atmospheric    pressure.  The 

dashed  curve  in  these  figures  represents  a  least-squares  fit 
for  data  taken  on  the  "optimum"  tube  with  rectangularly- 
shaped  fins  and  a  fin  height  of  1.0  mm  (tube  number  6).  In 
both    cases,      the    "Wolverine"    tube    (tube    number    43)       gave   the 
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worst  performance.  This  was  due  primarily  to  the  small  fin 
spacing  (0.7  mm)  that  resulted  in  a  completely  flooded  tube. 
Based  on  the  results  of  Georgiadis  [7],  the  small  fin  thick- 
ness (0.3  mm)  of  the  Wolverine  tube  may  also  be  responsible 
for  its  poor  performance.  Georgiadis  showed  that  for  a  fin 
spacing  of  1.0  mm  the  tubes  with  a  fin  thickness  of  0.5  mm 
(tube  number  16)  performed  16  percent  and  6  percent  poorer, 
under  vacuum  and  at  atmospheric  pressure  respectively,  than 
the  tube  with  a  fin  thickness  of  1.0  mm  (tube  number  5)  . 
Eased  on  this  observation,  the  small  fin  thickness  is  also 
responsible  for  the  poorer  performance  of  the  "Wolverine" 
tube  in  comparison,  for  example,  to  the  tube  with  fins  of 
parabolic  shape  (tube  number  37),  which  was  also  completely 
flooded  by  condensate.  Therefore,  the  poor  performance  of 
this  tube  is  probably  due  more  to  increased  condensate 
flooding  and  smaller  fin  thickness  than  to  fin  shape.  On 
the  other  hand,  the  best  performance  was  observed  for  the 
spirally-finned  tube  (tube  number  36)  under  vacuum  condi- 
tions, and  the  tube  with  parabolic  fins,  where  fin  spacing 
was  0.7  mm  (tube  number  37),  at  atmospheric  pressure.  The 
fact  that  tube  number  37  performed  well  at  atmospheric  pres- 
sure, but  poorer  than  the  spirally-finned  tube  under  vacuum 
conditions  was  somewhat  unexpected.  It  must  be  kept  in  mind 
that  these  data  runs  were  repeated  at  least  three  times  on 
different  days,  thus  confirming  this  trend.  One  possible 
explanation  centers  on  variable  fluid  properties,  as 
discussed  below. 

Under  vacuum  conditions,  the  condensate  retention 
angle  is  significantly  smaller  for  tube  number  36  (about  30 
degrees)  than  for  tube  number  37  (completely  flooded) ,  and 
the  saturation  temperature  is  about  50  K  lower.  Therefore, 
the  relatively  high  heat-transfer  coefficient  of  the 
spirally-finned  tube  is  due  mainly  to  its  small  condensate 
retention  angle.   When  the  condensing  temperature  is  raised, 
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both  surface  tension  and  viscosity  decrease.  For  example, 
the  surface  tension  decreases  by  about  10  percent,  and  the 
viscosity  decreases  by  about  45  percent  when  the  temperature 
is  raised  from  40  °C  to  90  °C  (assumed  approximate  film 
temperatures).  Decreased  viscosity  should  result  in  higher 
heat-transfer  coefficients  at  atmospheric  pressure  than 
under  vacuum  conditions  for  a  given  tube.  Further,  it  may 
be  possible  that  the  effect  of  viscosity  is  dependent  on  the 
fin  shape.  Unfortunately,  no  equations  exist  to  support 
this  hypothesis.  The  effect  of  surface  tension,  on  the 
other  hand,  is  fairly  well  understood;  the  condensate  reten- 
tion angle  decreases  with  decreasing  surface  tension. 
Therefore,  the  reduction  in  retention  angle  (from  vacuum  to 
atmospheric  pressure)  for  the  tube  with  "parabolic"  fins  is 
greater  than  the  reduction  for  the  spirally-finned  tube, 
since  the  latter  was  already  small  under  vacuum  conditions. 
Thus,  the  reduced  retention  angle  on  the  tube  with  "para- 
bolic" fins  provides  an  increased  heat-transfer  performance, 
surpassing  the  performance  of  the  spirally-finned  tube. 

The  performance  of  the  remaining  two  tubes  (tubes 
number  34  and  35)  is  the  same  under  vacuum  conditions,  but 
at  atmospheric  pressure,  the  tube  with  triangularly-shaped 
fins  (tube  number  34)  outperformed  the  tube  with  trapezoidal 
fins  (tube  number  3  5)  by  about  10  percent.  An  argument 
similar  to  the  one  above  may  explain  this  behavior. 

2.   Effects  on  Enhancement  Ratio 

Referring  to  Table  II,  it  can  be  seen  that  all  of 
the  above  tubes  provided  significant  heat-transfer  enhance- 
ments over  smooth-tube  values.  Note  that  the  enhancement 
obtained  under  vacuum  conditions  for  tube  number  36  was  the 
same  as  obtained  for  the  "optimum,"  rectangularly-shaped 
finned  tube  with  the  same  fin  height  (tube  number  6) .  Under 
atmospheric   conditions,   however,    the  tube   with  fins   of 
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parabolic  shape  (tube  number  37)  performed  significantly 
better  than  any  of  the  tubes  with  rectangular  fins  of  the 
same  height  (1.0  mm)  . 

like  the  tubes  with  rectangular  fins,  most  of  these 
tubes  displayed  enhancements  above  the  area  ratio,  as  shown 
by  the  two  Eo/Ar  columns  in  Table  II.  These  ratios  range 
from  0.94  to  2.95  under  vacuum  conditions,  and  from  1.52  to 
4. 13  at  atmospheric  pressure.  Like  tube  number  15,  the 
small  values  of  0.94  and  1.52  for  the  "Wolverine"  tube 
resulted  from  the  large  amount  of  retained  condensate,  and 
by  the  very  thin  fins  (about  0.3  mm)  as  discussed  earlier. 
As  was  the  case  for  rectangular  fins,  the  highest  Eo/Ar  of 
2.95  and  4.13  were  shown  for  the  tube  (tube  number  36)  with 
a  relatively  small  area  ratio  of  1.37,  indicating  once  again 
the  importance  of  surface-tension  effects. 

3 •   The  Performance  of  " Parabolic"  Fins 

In  order  to  study  the  parabolic  fin  more  closely,  a 
tube  with  parabolic  fins  was  manufactured  (tube  number  38) 
with  the  same  fin  height,  spacing,  and  thickness  as  tube 
number  17  (tested  by  Georgiadis  [7]  and  this  investigator). 
The  performance  of  these  two  tubes  is  compared  in  Figure 
5.12  under  vacuum  conditions,  and  in  Figure  5.  13  at  atmos- 
pheric pressure.  Clearly,  the  parabolic  fins  provide  a 
greater  enhancement  than  the  rectangular  fins.  The  reason 
for  this,  as  mentioned  earlier,  is  probably  the  continuous 
change  of  radius  of  curvature  (increasing  from  the  trough 
area  to  fin  tip)  for  the  parabolic  shape.  The  resulting 
pressure  gradient  along  the  fin  height  is  more  favorable 
(than  on  rectangular  fins),  which,  in  turn,  leads  to 
improved  drainage  from  the  fin  surface.  This  effect  was 
seen  most  clearly  on  tube  number  37.  Another  observation 
(that  contributes  to  improved  heat  transfer)  made  on  the 
tube   with  "parabolic"   fins  is   that,   although   completely 
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flooded,  the  fin  tips  (about  half  the  fin  height)  were 
visible  along  the  majority  of  the  fin  circumference  (except 
for  about  the  bottom  30  degrees) . 

As  discussed  in  section  E.1,  tube  number  37 
performed  poorer  than  the  spirally  finned  tube  (tube  number 
36)  under  vacuum  conditions.  The  reasons  for  this  are  not 
fully  understood;  so,  further  tests  of  these  fin  geometries 
(i.e.,  spiral  and  parabolic)  should  be  performed. 

P.   EFFECT  OF   FIN-METAL  THERMAL  CONDUCTIVITY   ON  STEAM-SIDE 
COEFFICIENT 

To  investigate  the  effect  of  fin-metal  thermal  conduc- 
tivity on  the  heat- transfer  performance,  a  stainless  steel 
tube  (tube  number  40)  was  manufactured  with  the  same  fin 
dimensions  as  the  "optimum"  copper  tube  (tube  number  6)  with 
rectangular  fins,  and  a  fin  height  of  1.0  mm.  Due  to  the 
low  thermal  conductivity  of  stainless  steel,  it  was  neces- 
sary to  reduce  the  tube  wall  thickness  to  0.5  mm  to  prevent 
an  excessive  wall  thermal  resistance.  A  copper  tube  (tube 
number  39)  with  the  same  dimensions  as  the  stainless  steel 
tube  was  manufactured  for  comparison.  Also,  two  smooth 
tubes  having  an  outside  diameter  equal  to  the  root  diameter 
of  the  finned  tubes  (tubes  number  41  and  42)  were  manufac- 
tured; the  first  tube  was  made  of  copper,  and  the  other  was 
made  of  stainless  steel.  The  results  for  data  runs  taker, 
under  vacuum  conditions,  and  at  atmospheric  pressure  are 
shown  by  Figures  5.14  and  5.15  respectively. 

As  expected,  the  copper  finned  tube  provided  a  signifi- 
cant enhancement,  while  the  stainless  steel  finned  tube 
actually  reduced  the  heat-transfer  performance  slightly  ir. 
comparison  to  the  smooth  tube  (tube  number  42) .  The  thermal 
conductivity  of  stainless  steel  is  much  lower  than  that  of 
copper   (i.e.,   15  W/m.K   compared   to   385  W/m.K),   so   the 
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thermal  resistance  of  the  stainless  steel  fins  is  signifi- 
cantly higher,  resulting  in  a  poorer  fin  efficiency-  This 
is  combined  with  the  deleterious  effects  of  retained  conden- 
sate (introduced  by  the  use  of  fins)  so  that  the  above 
factors  dominate  over  the  beneficial  thinning  of  the  conden- 
sate film  between  fins.  The  net  effect,  then,  is  a  slight 
degradation  of  the  heat- transfer  performance  compared  to  a 
smooth  tube. 

G.   EFFECT  OF  STEAM  VELOCITY 

As  discussed  in  Chapter  II,  most  models  of  film  ccnden- 
sation  on  finned  tubes  have  neglected  the  effect  of  vapor 
shear.  For  this  reason,  a  decision  was  made  to  perform  runs 
for  the  "optimum"  tube  (tube  number  6)  with  rectangular 
fins,  and  fin  height  of  1.0  mm,  at  nominal  steam  velocities 
of  2,  4,  6,  and  8  m/s.  Steam  velocity  through  the  test 
section  was  controlled  by  the  boiler  input  power,  and  since 
a  steam  velocity  of  1  m/s  was  the  maximum  obtainable  at 
atmospheric  pressure,  the  data  runs  were  performed  under 
vacuum  conditions  only.  The  results  of  data  taken  for  the 
above-mentioned  steam  velocities  are  shown  in  Figure  5.16. 

The  trend  shown  in  Figure  5.16  indicates  that  increased 
steam  velocity  (i.e.,  increased  vapor  shear)  leads  to 
enhancements  of  up  to  10  percent  as  steam  velocity  increases 
from  2  to  8  i/s.  This  is  to  be  expected  since  the  shear 
stress  at  the  liquid-vapor  interface  tends  to  pull  the 
condensate  toward  the  bottom  of  the  tube,  thereby  thinning 
the  condensate  film  at  the  top  part  of  the  tube. 

The  above-mentioned  data  were  compared  to  data  taken  by 
Yau  et  al.  [1]  for  steam  velocities  of  0.5,  0.7,  and 
1.1  m/s  at  atmospheric  pressure.  They  reported  enhancements 
of  up  to  40  percent  for  the  steam  velocity  of  1.1  m/s 
compared   to   the   0.5  m/s    velocity.     The   cather   large 
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difference  in  enhancements  reported  by  Yau  et  al.  compared 
to  those  of  this  investigation  was  not  expected.  It  should 
he  kept  in  mind,  however,  that  the  data  of  Yau  et  al.  were 
collected  at  atmospheric  pressure  (compared  to  vacuum  condi- 
tions for  this  investigation),  where  the  increased  satura- 
tion pressure  leads  to  a  lower  condensate  viscosity  (about 
45  percent  as  discussed  in  section  E. 1) .  This,  in  turn, 
allows  condensate  to  flow  more  freely  so  that  the  effect  of 
shear  stresses  acting  on  the  liquid-vapor  interface  would  be 
more  pronounced.  Another  factor  to  be  considered  is  the 
slightly  larger  fin  spacing  used  by  Yau  et  al.  (2.0  mm 
compared  to  1.5  mm  for  this  investigation) .  If  the  fin 
spacing  is  very  small,  the  vapor  shear  may  not  be  "felt"  on 
the  sides  of  tha  fins  or  on  the  fin  bases.  This  is  due  to 
the  fact  that  the  steam  flows  around  the  tube,  thus  leaving 
nearly  stagnant  steam  between  fins.  As  the  fin  spacing 
increases  though,  the  fin  sides  and  bases  will  be  more 
"accessible"  to  steam,  therefore,  giving  rise  to  the  benefi- 
cial effects  of  vapor  shear.  Thus,  the  larger  enhancements 
of  Yau  et  al.  [ 1  ]  compared  to  this  investigation  may  be 
explained  by  the  two  phenomenon  just  described. 

An   attempt   was  made   to   correlate   the  data   using   a 
Fujii-type  [29]  equation  given  by: 


Nu 


=  c  F 


;H 


Rennu'b  (5.3) 


r-P 


where 


Hu  =  haf  Do  (5.4) 
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Re 


tP 


=    Pf    Vs    no 
Mf 


(5.5) 


F   =      g   Do    Mf   hfR 
Vsz   kf    AT 


(5.6) 


V      =    steam    velocity 


"af 

EetP 

c,    d 


steam-side    heat-transfer    coefficient   based 
on    actual   surface   area    (including    fins) 

two-phase   Reynolds    number 

constants    to    te    determined    experimentally 


The  constants  c  and  d  used  in  the  Fujii  equation  [29],  which 
was  developed  to  account  for  vapor  shear  on  smooth  tubes, 
were    0.96      and    0.2    respectively.  Values   of      c   =    1.30  2   and 

d  =  0.236,  however,  were  found  to  more  closely  fit  the  data 
of  this  investigation  after  a  least-squares  fit  was 
performed.  A  comparison  of  the  experimental  data  to  this 
curve    is      shown    in      Figure    5.17.  Although    a      fairly    close 

agreement  has  been  obtained,  the  orientation  of  data  clus- 
ters does  not  match  the  slope  of  the  least-squares  line 
generated  by  globally  fitting  all  the  data.  It  is  possible 
that  the  errors  in  the  experimental  constants  used  in  the 
Sieder-Tate-type  equation  are  partly  responsible  for  this 
questionable  trend.  Georgiadis  [7]  showed  that  using  larger 
values       (from         the      "modified        Wilson      method")  of      the 

Sieder-Tate  coefficient  (C)  tends  to  rotate  the  data  clus- 
ters to  more  closely  match  the  experimental  correlation. 
Further,  a  Fujii-type  correlation  may  be  inappropriate  for 
finned   tubes    due    mainly    to   condensate    retention    that    occurs. 
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Note  that  the  slope  of  the  line  predicted  by  equation 
(5.3)  for  the  experimental  data  does  not  match  the  slope  of 
the  Fujii  equation  mentioned  earlier,  but  instead,  nearly 
matches  the  slope  of  the  Nusselt  equation  for  smooth  tubes 
given  by: 


h   =  0.728 


3       2  ,1/4 

kf     Pf    fl    hft> 
|if   n0    AT 


(5.7) 


which    can   be    rewritten   as 


Nu 


!/«• 


FT [TTT 

Ketn 


=    0.7?R    F 


(5.8) 


In  this  manner,  the  effect  of  vapor  shear,  at  least  on  the 
tube  tested  (tube  number  6) ,  is  very  small  compared  to 
smooth  tubes  [7].  In  fact,  using  the  experimentally-found 
d-value  in  equation  (5.3),  the  steam-side  coefficient  and 
steam    velocity    are    related    as    follows: 


h  ~  v 


0.028 


(5.9) 


Once  again,  this  observation  shows  that  the  effect  of  steam 
velocity  on  the  steam-side  heat-transfer  coefficient  is  very 
small. 

A  vertical  shift  of  the  least-squares-fit  compared  to 
the  Nusselt  theory  results  from  the  insreased  c-value  (1.302 
compared  to  0.728)  of  equation  (5.8) .  In  fact,  the  ratio  of 
the  experimental  c-value  to  the  coefficient  in  the  Nusselt 
equation  (i.e.,  1.302/0.728  =  1.788)  represents  the  heat- 
transfer    enhancement    in      excess   of    the    area      increase   due   to 
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finning.      Therefore,      the   c-value   depends   fin  geometry,      and 
surface   tension    of    the  condensate   etc.. 

In  view  of  the  above  discussion,  further  attempts  to 
modify  the  Nusselt  eguation  to  include  vapor  shear,  surface- 
tension  effects,  and  fin  geometry  etc.  may  lead  to  an  accu- 
rate correlation  for  the  steam-side  heat-transfer 
coefficient    for    finned   tubes. 

H.       EFFECT      OF    INTERNAL      AND      EXTERNAL       ENHANCEMENTS    ON      THE 
OVERALL    COEFFICIENT 

Up  to  this  point,  all  discussions  have  delt  with 
enhancements  to  the  steam-side  heat-transfer  coefficient. 
Even  with  large  steam-side  enhancements,  however,  a  large 
thermal  resistance  through  the  tube  wall,  or  on  the  water 
side  may  dominate,  and  completely  negate  the  steam-side 
enhancement.  For  this  reason,  a  decision  was  made  to  test 
the  "optimum"  tube  with  rectangular  fins  (fin  height  of 
1.0  mm  -  tube  number  6),  and  the  corresponding  smooth  tube 
(tube  number  1) ,  with  and  without  the  insert.  The  results 
of  these  data  runs,  taken  under  vacuum  conditions,  are  shown 
in  Figure  5.18.  At  a  water  velocity  of  about  3  m/s,  use  of 
the  insert  in  the  smooth  tube  resulted  in  a  40  percent 
increase  of  the  overall  coefficient,  as  did  the  use  of 
external  fins   without  the   insert.  By   combining    the    use  of 

fins  and  the  insert,  an  enhancement  of  nearly  155  percent 
was  attained.  It  must  be  noted,  however,  that  the  insert 
used  during  these  experiments  is  not  suited  for  use  in 
condensers  due  to  the  very  high  pressure  drop  it  creates. 
Instead,  Webb  et  al.  [30]  showed  that  the  use  of  internal 
flutes,  fins,  or  circumferential  ribs  may  be  considered  to 
increase  the  water-side  coefficient  by  a  factor  of  up  to  2.0 
(compared  to  2. 5  obtained  by  the  insert  used  during  this 
investigation)  .         As   a      result,      it    appears   that      the   use  of 
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external  fins  togather  with  internal  enhancements,  could 
increase  the  overall  heat- transfer  coefficient  by  approxi- 
mately   100    percent. 

I.       DEVELOPMENT    OF    AH   EXPERIMENTAL    CORRELATION 

Using  the  1985  model  of  Webb  et  al.  [25]  (discussed  in 
Chapter  II)  as  a  starting  point,  attempts  were  made  to 
correlate  the  experimental  data  obtained  for  tubes  with 
rectangularly-shaped    fins.  Time    limitations      prevented   an 

in-depth  effort,  but  a  few  simple  modifications  were  made. 
Before  discussing  these  modifications,  the  general  scheme  of 
the  Febb  et  al.  model  will  be  discussed.  A  more  detailed 
step-by-step  discussion   is   provided    in   Appendix    A. 

First,  a  value  for  the  parameter  £  must  be  determined 
for  the  fin  geometry  in  question.  With  this  value  known, 
the  heat-transfer  coefficient  for  the  unflooded  fin  surfaces 
can  be  obtained  using  equation  (2.30).  Values  for  the  heat- 
transfer  coefficients  of  the  unflooded  tube  area  between 
fins,  and  for  the  flooded  region  are  determined  next  using 
equations    (2.32)         and      (2.33)       respectively.  The      above- 

mentioned  heat- transfer  coefficients  are  then  substituted 
into  equation  (2.31)  to  obtain  the  average  heat- transfer 
coefficient    for    the    finned    tube   in    question. 

A  careful  review  of  the  Webb  et  al.  model  revealed  two 
deficiencies  worth  noting.  First,  in  their  two-dimensional 
analysis  of  the  flooded  tube  region,  the  heat-transfer  coef- 
ficient for  the  flooded  root  area  between  fins  was  applied 
to  the  fin  tips.  Figure  2.U  shows  that  a  constant,  rela- 
tively thick  condensate  film  is  assumed  which,  when  applied 
to  the  fin  tips,  results  in  an  underprediction  of  the  heat 
transfer.  Ose  of  the  Nusselt  equation  for  smooth  tubes 
would  be  more  appropriate  for  fin  tips  in  this  region. 
Additionally,       the    effects    of    convection      in    this    region   are 


95 


& 


"—    u 
<D    to 


w 

tn 

\ 

X 

a 

e 

s 

o 

_a                                     • 

m 

CN 

CO 

II 

ii 

U] 

0) 

> 

a, 
1 

<D 

-P 

•P 

XI 

M 

M 

3 

0) 

0) 

E-< 

U3 

0) 

en 

C 

XI 

c 

£ 

H 

a 

M 

+J 

Eh 

0 

A 

■U 

0 

P 

JS 

3 

in 


OJ 


CD 


ID 


CO 


/•s 

(0 

\ 

E 

> 


N 


0) 


(D 


CO 


O 


O 


u 

O 

«H 

4-) 

a 

0) 

•H 

U-fl 

•H4J 

4-t-H    • 

W  *  — 

(V      (0 

one 

UX5  3 

90S 

UH 

a       • 

4-113  0 

W  0)  (d 

(3  et> 

(TJ  Cw 

H-H 

H&h+J 

1        M 

•P^  <D 

to  a  w 

0>  n3  C 

DC      M 

a) 

HJ-W 

HAS 

ASH  O 

M      A 

(UJ54J 

>4->.H 

O  O* 

0 

flrO 

in  c 

oo     (0 

*~  (ti 

• 

m 

0) 

u 

=j 

D» 

•H 

fM 

c*  ■a.w/wH)    °n 


96 


ignored;   however,    to  include  these   effects  would   make  a 
solution  extremely  complex. 

The  second  deficiency,  probably  valid  only  for 
rectangularly-shaped  fins,  involves  the  assumed  condensate 
film  shape  between  fins  on  the  unflooded  tube  surface. 
Figure  2.4  shows  a  depression  in  the  film  next  to  the  fin 
base,  with  an  increasing  condensate  thickness  away  from  it. 
Condensate  film  shapes  observed  during  this  investigation 
continued  to  decrease  away  from  the  fin  base  as  shown  by 
Figure  2.1,  probably  due  to  the  thinning  effect  of  surface- 
tension  forces.  Based  on  the  assumption  of  a  relatively 
thick  condensate  film,  predictions  using  the  Webb  et  al. 
model  will  underpredict  heat  transfer  from  this  part  of  the 
tube. 

To  use  the  Webb  et  al.  model  for  rectangular  fins,  a 
value  of  -0.9  was  chosen  for  e,  .  A  value  of  -1.0,  which 
corresponds  to  a  straight  line  (i.e.,  straight  side  of  the 
rectangular  fin)  would  seem  more  appropriate;  however,  this 
value  would  lead  to  zero  heat  transfer  from  the  fin  sides 
using  equation  (A.1).  Webb  et  al.  used  values  around  -0.7 
to  -0.8  for  their  slightly  trapezoidal  fins;  so  the  value  of 
-0.9  was  assumed  to  be  reasonable.  Using  this  value  along 
with  a  <f>  value  of  5.0  (used  by  Webb  et  al.)  resulted  in 
significant  underprediction  of  the  experimental  data 
obtained  by  this  investigator  and  by  Georgiadis  [7],  as 
shown  in  Figure  5.19. 

In  an  attempt  to  improve  the  predictions,  values  for  the 
<j>  ratios  were  calculated  using  the  paralled  heat  flow  paths 
proposed  by  Owen  et  al.  [13]  discussed  earlier  ir. 
Chapter  II.  These  4>  values  are  dependent  on  fin  spacing, 
and  should  decrease  with  increasing  fin  spacing.  This  trend 
was  observed,  and  for  the  set  of  tubes  with  both  fin  height 
and  fin  thickness  of  1.0  mm  (s  =  0.5,  1.0,  1.5,  2.0,  and 
U.O  mm  -   tubes  number  1  to   8),   resulted  in  4>   values  that 
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ranged   from    1.7     (s    =    0.5    mm)     to    1.2     (s    =    U.O    ram)  .       Using    the 
Febb    et    al.       model    once    again,    tut    with    a    constant    average 
value   of    1.5,       gave    the    results   shown    in    Figure    5.20    for    the 
above-mentioned    tubes    (tubes      number    4    to    8) .  This    Figure 

also  shows  a  plot  of  the  Webb  et  al.  model  using  the  vari- 
able (i.e.,  1.7  to  1.2)  which  made  almost  no  difference. 
Cnce  again,  the  experimental  data,  also  shown  in  Figure  5.20 
were    significantly    underpredicted. 

In  order  to  correct  for  the  assumed  thick  condensate 
layer  between  fins,  the  work,  of  Fujii  et  al.  [31]  was  used. 
Their  paper  delt  with  smooth  tube  enhancements  from  coiling 
a  wire  around  the  tube  circumference,  and  provided  an  equa- 
tion to  account  for  the  thinning  effect  between  wire  wraps 
as    shown  below: 


Mu 


f  . 


Mu, 


(s  +  t) 


4    (I  +  A) 


,    1/- 


(5.  10) 


where 


A   =         j    °  Do 

~ —T~ 


(5.11) 


Pf    R    s 


r      =   0.03 
s 


2    a 


pf  r. 


3/2 


-2 


w 


(5.12) 


D      =    wire   diameter    (taken   as    the    fin    height 


w 


for  use  in  the  Webb  et  al.  model) 


Adding  this  correction  factor  to  the  heat-transfer  coeffi- 
cient for  the  unflooded  root  area  between  fins  made  a 
significant   improvement  for   the  group   of  tubes   mentioned 
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earlier  (tubes   number  4   to  8)    as  shown   by  Figure   5.20. 
Nonetheless,  experimental  data  were  still  underpredicted. 

In  summary,  it  should  be  kept  in  mind  that  this  model 
was  not  developed  for  rectangular  fins.  The  most  signifi- 
cant discrepancy,  however,  appears  to  be  on  the  prediction 
of  heat  transfer  through  the  flooded  area.  This  model 
predicts  only  minimal  heat  transfer  through  the  flooded 
portion  of  the  tube,  whereas  Fanniarachchi  et  ai .  [2] 
reported  enhancements  at  least  equal  to  the  area  ratio  for 
completely  flooded  tubes.  Additionally,  there  may  be  other 
phenomena  which  have  not  been  taken  into  account  (such  as 
intense  convection  in  the  condensate)  which,  if  included, 
may  improve  the  predictions.  Nonetheless,  this  model  shows 
promise  and  should  be  pursued  further. 
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71.    CONCLUSIONS    AND    RECOMMENDATIONS 


A.       COHCLOSIONS 


The  use  of  fins  may  lead  to  significant  enhancements 
of  the  steam-side  heat -transfer  coefficient.  Under 
vacuum  conditions,  enhancements  as  high  as  4.0  (tube 
number  36)  were  realized,  while  enhancements  of  up  to 
6.9  (tube  number  37)  were  realized  at  atmospheric 
pressure.  These  enhancements,  which  were  up  to  three 
times  greater  than  the  area  ratio  (finned  tube 
area  /  smooth  tube  area)  are  due  to  the  thinning 
effect   of   surface-tension    forces. 

For  rectangularly-shaped  fins  with  fin  heights  of  0.5 
and  1.5  mm,  the  2.0  mm  fin  spacing  appears  to  be  the 
optimum  among  the  fin  spacings  tested  (1.0,  1.5,  2.0, 
and  4.0  mm).  This  was  the  case  for  both  vacuum  and 
atmospheric  conditions.  For  similar  tubes  with  fin 
heights  of  1.0  and  2.0  mm,  Georgiadis  [7]  found  the 
optimum   spacing   to   be    1.5    mm. 

The  use  of  fins  on  a  stainless  steel  (low-thercal- 
conductivity- metal)  tube  degraded  the  heat-transfer 
performance  in  comparison  to  its  corresponding  smooth 
tube.  In  this  case,  the  poor  fin  effeciency  combined 
with  the  deleterious  effects  of  retained  condensate 
dominated  over  the  beneficial  thinning  of  the  conden- 
sate film  between  fins,  resulting  in  a  net  reduction 
of  the   heat- transfer   performance. 

Fin  geometry  affects  the  extent  to  which  the  surface- 
tension  forces  thin  the  condensate  film  on  the  fin 
surfaces.         Parabolically-shaped      fins   were      shown    to 


102 


maximize   the   thinning    effect   of   surface-tension 
forces. 

5.  Increasing  steam  velocity  (i.e.,  vapor  shear)  from  2 
to  8  m/s  increased  the  steam-side  coefficient  by 
about  10  percent  (under  vacuum  conditions) ,  in 
contrast  to  Yau  et  al.  [ 1 ],  who  showed  a  40  percent 
increase  when  steam  velocity  was  increased  from  0.5 
to  1.1  m/s  (at  atmospheric  pressure). 

6.  For  tha  "optimum"  tube  with  rectangular  fins  and  fin 
height  of  1.0  mm  (tube  numbar  6)  ,  the  overall  heat- 
transfer  coefficient  without  the  insert  showed  an 
enhancemant  of  40  percent  over  the  corresponding 
smooth  tube.  With  the  insert  in  place,  the  finned 
tube  showed  overall  enhancements  of  about  155 
percent. 

7.  Most  theoretical  models  developed  thus  far  to  predict 
the  heat-transfer  performance  of  steam  condensation 
on  horizontal  finned  tubes  neglect  one  or  more 
phenomena  of  major  importance.  Thus,  there  are  no 
models  currently  available  to  accurately  predict  the 
performance  of  finned  tubes  used  to  condense  steam. 

B.   RECOMMENDATIONS 

1.  Make  further  refinements  to  the  Webb  et  al.  [25] 
model  in  an  effort  to  predict  the  experimental  data 
within  a  reasonable  accuracy. 

2.  Manufacture  and  test  additional  finned  tubes  with 
parabolic  and  spiral  fin  shapes  with  different  dimen- 
sions (i.e.,  fin  spacing,  haight,  and  thickness) .  In 
this  manner,  determine  the  optimum  dimensions  for 
each  fin  shape. 
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3.  Manufacture  and  test  a  series  of  tubes  with  varying 
tube  diameter  to  study  the  effect  of  diameter  on  the 
heat-transfer  performance.  For  this  purpose,  select 
the  optimum,  rectangular-shaped  fins  (s  =  1-5  mm, 
e  =    1.0    mm,    and   t  =    1.0    mm). 

4.  Collect  data  on  finned  tubes  fitted  with  drainage 
strips  made  out  of  both  solid  and  porous  metal  to 
determine   the   enhancements   they   provide. 

5.  Continue  tests  on  finned  tubes  with  tube-metal 
thermal  conductivities  between  those  of  copper  and 
stainless  steel,  such  as  copper-nickel  alloys  and 
aluminum. 

6.  Conduct  data  runs  at  atmospheric  pressure  with  steam 
velocities  of  0.5,  0.7,  and  1.1  m/s,  using  tube 
number   6    to    compare  with  the    data   of   Yau   et   al.    [1]. 

7.  Perform  data  runs  using  different  fluids  with  widely 
varying  fluid  properties  (i.e.,  surface  tension, 
viscosity,  etc.)  to  further  study  these  effects,  and 
to  compare  with  existing  theoretical  models  and 
experimental    correlations. 
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APPENDIX    A 
PROCEDURE    FOR    DSE    OF    THE    WEBB    ET    AL.     [25]    MODEL 

As    proposed    by    Webb    et      al.,      solving    for    the    steam-side 
heat-trar.sfer    coefficient    proceedes    as    follows: 

1.  First,  a  value  for  £  must  be  determined.  This 
parameter,  as  shown  by  figure  2.4,  is  a  measure  of 
the  fin  aspect  ratio  (e/t) ,  and  it  must  be  deter- 
mined using  an  iterative  method.  It  is  a  function 
of  fin  geometry  only,  and  should  match  the  fin  side 
profile  as  closely  as  possible.  Webb  et  al.  [25] 
provide  details  of  this  method,  and  no  further 
discussion  is  provided  on  this  since  a  true  ^-value 
does  not  exist,  as  defined  by  Adamek  [27],  for  a 
rectangular-shaped    fin. 

2.  With  the  value  of  known,  determine  the  heat- 
transfer  coefficient  for  the  fin  sides  in  the 
unflooded  portion  of  the  tube.  This  value  is 
obtained    using    the    following    equation    [25,    27]: 


hf    =    2.140 


q   hfp    °m   Sm    Pf    fg  +    1) 

Mf  Iff  at  7TTTF1" 


L/H 


(A.1) 


Determine  the  heat-transfer  coefficient  for  the 
unflooded  root  area  between  fins  using  the  Nusselt 
equation  written  in  terms  of  the  film  Reynolds 
number.  For  this  purpose,  use  the  following  itera- 
tive  nethod: 
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a)       Calculate    the    condensation    rate     (m)     using: 


m   =   n   hf   Apr    AT 
hfp 


(A. 2) 


b)       Assume: 


m      =   m 

r 


(A. 3) 


c)       Calculate   the    film   Reynolds    number    by: 


Ref    = 


4    n, 


pf  rr.s  +  to  -  tb] 


(A. 4) 


d)       Calculate   hr    using 


h      =    1.514 

r 


2 
»*    ,       Re- 


kf  J    pf^    f» 


-1/3 


(A. 5) 


e)       Calculate    the    base   heat    flow   by: 


0     =   h     A,  ..    AT 
r  r      ht 


(A. 6) 


f)       Solve  for   m    using: 


m   =  m  + 
r       h 


fg 


(A. 7) 


g)   Repeat  steps  c  thru  f  until  convergence  is 
obtained. 
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Determination  of  the  heat-transfer  coefficient  in 
the  flooded  part  of  the  tube  requires  that  a 
computer  solution  be  used  to  solve  the  two- 
dimensional  conduction  problem  of  heat  flow  through 
the  fins  and  retained  condensate.  This  solution 
provides  a  ratio  (tf>)  of  actual  heat  flow  into  the 
coolant  divided  by  the  heat  flow  that  would  result 
from  fins  of  zero  thickness.  The  heat- transfer 
coefficient    is    then    determined    by: 


\  = 


]1L 


(A. 8) 


The      fin    efficiency      is      determined      next    using      the 
following   equation: 


n  = 


tanh    fn    (e   +   t  /2  ^  ] 
n    (e  +  t/2) 


(A. 9) 


where 


m   = 


hf 


1/2 


(A.  10) 


It  should  be  kept  in  mind  that  this  fin  efficiency 
is  only  valid  for  "thin"  fins.  Nonetheless,  its  use 
is  justified  since  there  are  no  equations  available 
specifically    for    "thick"    fins. 
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6.  The  heat-transfer  coefficient  for  the  entire  tube, 
then,  is  obtained  by  substituting  the  above  quanti- 
ties  into    the   following   equation: 


h       =  h  n     - 
ow  o 


Asf 


+    n   h 


Aft' 

TT     -      |p 

+     h,      ^ 

Asf  _ 

7T 

b  * 

(A.  11) 
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APPENDIX    B 
MODIFIED    HILSON    METHOD 

The      "modified       Wilson    method"      calculates      the      leading 
coefficient      for      the  Sieder-Tate      equation.  This      method 

assumes  a  form  of  correlation  for  both  the  water-side  and 
steam-side  heat-transfer  coefficients,  each  of  which 
contains  a  coefficient  to  be  determined  by  iteration.  In 
the  past,  a  Musselt-type  equation  was  used  for  the  steam- 
side    heat-transfer    coefficient    as    shown    below   [7]: 


h     = 
o 


kf    Pf    Cpf    -    pv^    hfj 


1/3 


(B.1) 


Equation  (B.1)  results  in  a  6  value  of  0.655  for  conditions 
where   no      vapor    shear      is    present.  The    presence      of    vapor 

shear  generally  results  in  a  higher  value,  which  must  be 
determined    iteratively.  The    use    of    this      equation   should, 

therefore,  be  restricted  to  relatively  low  vapor  shear 
conditions  (i.e.,  vapor  velocity  less  than  1.0  m/s)  .  Note 
that  a  Nusselt-type  equation  does  not  represent  the  "correct 
form,"  as  it  contains  no  terms  involving  vapor  shear.  To 
avoid  this  deficiency,  a  correlation  developed  by  Fujii  and 
Honda  [29]  which  takes  the  effect  of  vapor  shear  into 
account    was   used,    and   is   shown   below: 


Nu 
Re,.    u'b 


1/5 


=   0.06    F 


(B.2) 
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Rewriting  equation  (B.2)   to  express  h  as  a  function  of  heat 
flux  and  vapor  velocity  results  in: 


h      = 
o 


f?    hfs 


!/«♦ 


3       3        3 
"f   Do    pf   Vs 


m    1/8 


kf  =  b  r 


(B.3) 


The    constant      3       in    equation    (B.3)  is    to    be      determined   by 

iteration.  A    Sieuer-Tate- type      equation    is      used    for      the 

water-side   hea t- transfer    coefficient,    as    shown   below: 


_± L  =    C   Re  Pr 


V.t 


0.  14 


=  c  n 


(B.4) 


Equations  (B.3)  and  (B.4)  are  then  substituted  into  the 
equation  for  the  overall  heat-transfer  resistance  given 
below  : 


TL.    A 


o         N    A1         h0    An 


(B-5) 


A   linear   equation  used   to   generate   the  Wilson   plot   is 
obtained,  and  shown  below: 


Y  = 


*  +  I 

C   3 


(B.6) 


v  = 


(B.7) 


Y  - 


*+i 

C        3 


(B.8) 
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The  parameters  X  and  Y  are  determined  from  fluid-property 
values  and  heat-flux  measurements  taken  during  data  runs 
using  an  uninstrumented  smooth  tube.  Iterations  between  the 
coefficient  C,  and  the  coefficient  8  ,  are  continued  until 
convergence  of  the  coefficients  (between  two  successive 
iterations)  within  0.1  percent  occurs.  The  slope  of  the 
Wilson  plot  generated  is  the  reciprocal  of  the  desired 
Sieder-Tate  coefficient. 
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APPENDIX  C 
LISTING  OF  RAW  DATA 

The  following  pages  contain  raw   data  obtained  for  tabes 

number  26  thru  44  under  vacuum  conditions  and  at  atmospheric 

pressure.   Raw  data  for  tubes  number  1  thru  25  are  presented 
by  Georgiadis  [  7 ]. 
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APPENDIX    D 
UNCERTAINTY    ANALYSIS 

There  is  always  an  uncertainty  associated  with  any  meas- 
urement which  is  dependent  on  the  measuring-device  accuracy, 
as    well      as   on    the      operator's   experience.  Numerical   data 

collected  during  this  thesis  effort  were  used  together  with 
theoretical  formulations,  so  final  values  of  the  steam-side 
heat-transfer  coefficient  may  be  distorted  due  to  uncer- 
tainty propagation  during  calculations.  In  cases  where  the 
final  results  show  large  uncertainties,  it  may  be  unwise  to 
draw  any  conclusions.  Instead,  the  apparatus  and/or  the 
measuring  techniques   may   have   to    be    modified. 

Estimates  of  uncertainties  for  this  investigation  were 
obtained  using  program  "UNA6"  which  is  listed  at  the  end  of 
this  appendix,  along  with  selected  uncertainty  evaluations. 
This  program  determines  the  uncertainties  using  an  equation 
proposed   by   Kline   and   McClintok   [32]    shown    below: 


3R 
~d7, 


-,    1/2 


(D.1) 


where 


WR  is  the  uncertainty  of  the  desired  dependent 
variable 

x,  ,  x2 ,  ....,  xn  are  the  measured  (independent) 
variables 

v.  ,  w~,  ....,  w   are  the  uncertainties  in  the  measured 
12  n 

variables 
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A  complete  discussion  covering  the  developement  of  the 
uncertainty  analysis  used  for  this  investigation  is  given  by 
Georgiadis    [7]. 
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r.:'i' 
( i i  0 ! 

020! 
930* 
040 

05  U 
060 
C7'J 
080 
090 
100 
no 
120 
130 
MO 
150 
160 
'70 

w: 

130 
190 
!31 
192 
200 
2!0 
220 
230 
240 
250 
250 

230 


-~3LZ  IMml 

rev::;::  : 


;NA6 

"en  r mat  v  2i 


COM  /Cc'  2(7) 

DIM  E<4> 

DAiH  0  .  10086091 


7 . 94369 . -  767245 . 3295 . 73025595 . 8 ! 


DATA  -9247486539.6. 97638E+ 11 . -2 .66' 92E* 1 3 . 2 . 940  7SE+1 4 
:<EAD  C(*> 

PRINT  USING  "lUX. ""DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 

PRINT 

PRINTER  IS  701 

BEE? 

INPUT  "ENTER  FILE  NAME".File5 

PRINT  USING  "10X.""File  Name:  "" . 12A" ;File£ 

BEEP 

INPUT  "ENTER  DaTA  SET  NUMBER  FOR  UNCERTAINTY  ANALYSIS" . Ids 

BEEP 

INPUT  "ENTER  PRESSURE  CONDITION  iO-V.!«A>  ".Pre 

Prc»Drc+1 

ASSIGN  *File  TO  ~ -lei 

ENTER  *Fi ie: Ifg.Inn 

IF  Ifg-U  THEN  ENTER  *File:Dd 

IF  Ifg-1  THEN  ENTER  SFi le :Dd .Dd.Dd 

FOR  1-1  TO  las 

ENTER    i»Fiie:Bvol  .  3amp  . Vtran.£tp.E<->  . rm.  7c  i  .  Tco.Phg.Pwater 

NEXT    I 

Ent =E<  9  > 

I!7    Prc=1     THEN 

3EEP 

PRINT    USING    "MUX  .""°ressure    Condition: 

Else 

PRINT  USING  "'OX .""Pressure  Condition. 


Vacuum    <  1  1    if  Pa) 

Atno^pneric    ('01    kPa)' 


2?u 

END  IF 

300 

PPINTER  IS  ' 

:;o 

BErD 

320 

^RINT  USING  "4X." 

'"Select 

tube  .i»a  i  :  ryoe  :  0 

330 

INPUT  Itt 

340 

IF  Itt-U  THEN 

350 

Do-. 01905 

560 

Di- .0127 

"".  "7  f"i 

r-l     s^r- 

c70 

tLJC. 

280 

Do-. 0145 

29o 

Di-  .0135 

400 

END  IF 

410 

PRINTER  IS  1 

420 

SEEP 

43U 

PRINT  USING  "ax . ' 

'"be  lee ' 

material  code: 

440 

PRINT  USING  "4X.' 

'"0-Cii. 

1  -  JO 

^50 

INPUT  inc 

460 

IF  lmc=0  THEM 

4  70 

K'c-335 

471 

Dfcc-10 

4f!0 

Ci-,071 

4^0 

Ac«0 

500 

ELSE 

510 

Kc-15 

511 

Dkc-I 

520 

C1-.O688 

530 

Ac  =  0 

540 

END  IF 

'tnm' 
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J. 


!c,50  fRINT£3  [3  70! 

1560  Ts-FNTvsv<Enf > 

[570  PRINT  USING  " I  OX . ""Steam  Temperature            -  "" . 3D . _D  .  ""  <Deg  C>' 

T-3 

1530  PRINT  USING  "'OX. ""Water  Fiou  Pace  U>           •  "". 3D. 2D":Fn 

1590  Dtci-.on: 

1600  Dtco-.005 

1610  BEEP 

1620  Dent -I  .OE-6 

T  630  Dts>SGR<((C<!  >+2»C(2>*Emf  *3*C<  3>-£nf  '2+4»0<4>-»Eraf  '  3>*D<?nf  >  2) 

16^0  T»< Tc:+Teo>/2   !  FILM  TEMPERATURE 

1650!  UNCERTAINTY  [N  THE  COOLING  WATER 

1660  Drho-.5  !  ERROR  IN  WATER  DENSITY 

1670  Dmt-.O08  !  ERROR  IN  MPSS  F-.QH  RATE 

1630  Rho-FMRho(T)  !  WATER  DENSITY 

1630  Mf-1  .048U5E-2+6.30932E->trT,   •  NASS  FLOl-1  RATE  OF  COOLING  WATER 

1700!  CORPECT  MF  FOP  [HE  TEMPERATURE  EFFECT 

1711]  (if  =  Mt"  *<  !  .0365-  i  .%64uE-3*'c  L»s.  252E-o*rci    2'/   39543a 

•  720  Ai*(PI*Di  2>/4  !  TUBE  INSIDE  CROSS  SECTION  AREA 

1730  Ddi-. 000025 

1740  "DaiJ<PI-0di  2>/4  »  ERROR  OF  INSIDE  TUBE  CROSS  AREA 

1750!  COMPUTE  THE  WATER  VELOCITY 

!760  v.j-Hf  /'Rho-Ai  >  !  WATER  VELOCITY 

'770  PRINT  USING  "I  OX  .""'Water  Yeiocitv              ■   "".Z.DD.""  fm/s) ;V 

1700!  CORRECT  OUTLET  *ATER  "E'-P .  rCR  THE*  FIXING  CHAMBER  EFrEC 

1700  IF  Inn-1  OR  Inn  =  5  r-IEM  Tco- Tco- .00^«Wy  C 

1600  IF  Inn-0  'HEN  "co-7co-<  -  .00  I  33*-  .001  »W«y  2) 

1  SI  0  T«<Tci*Tco>».5  '  FILM  TE*PF=?CTURE 

1S20!  COMPUTE  THE  ERROR  IN  *ATER  VELOCITY 

1320  Dvu»lv/.j«SQR<  (Dmt /Mf  >  2  +  <  Drno/Rho  >  2+<Dai/Ai>  2) 

IS40!  UNCERTAINTY  IN  THE  REYNOLDS  NUMBER 

1350  r1wFNMu<T>  !  WATER  VISCOSITY 

1060  Dmw-A.E-6  !  ERROR  OF  WATCR  VISCOSITY 

1370  Rp-<Rho»V.a*Di  >  /Mi.j 

1330  Dre-Re-SQR(<Drho/Rho> '2  +  <D./fi»/'Aj)  2*<Ddi/Di)  2+<Dou'Mu)  2> 

1390!  UNCERTAINTY  IN  THE  HEAT  TRANSFERRED 

1 300  Cpii»-FNCpu<T) 

1910  Q=Mf»<  Tco-Tci>«Cpiu 

1920  Dcpu-6 

'330  Dq  =  Q-SGR((Dnt"/Mf  ) -2^<  <Dtco/ <  Tcc-Tci)  >  >    2*-<  <  Otc:/<  Tco-Tc  i >  > )    2  +  <Dcp'j/Cpu>    2 

) 

1940!  UNCERTAINTY  IN  rHE    -EAT  -LUX 

1950  Dl-  .0005  !  ERROR  IN  tijbe  LENGKT 

I960  Ddo-. 000025 

1970  L-. 13335  !  CONDENSING  TURE  LENGTH 

1930  Qp-Q/(PI-Do-!_)  !  HEAT  FLUX 

1390  PRINT  USING  "'OX. ""Heat  -lux                   -  "".Z.2DE.""  <W/m  2) : 

Qp 

2000  Dap-<Qp/PI>*SQR<<Dq/Q)  2*<  Dao/Do>  '2*<  D  I/L>  2> 

2010  Lmtd=<  Tco-Tci)/L0G(  ( Ts-T.:i  >/<  T«s-Tco> > 

2020  Uo-Op/Lmtd  !  QVERALu.  HEAT  TRANSFER  COEF. 

2030  A1»0ts»< Tci-Tco)/< ( Ts- :ci)-< Ts-Tco>*LOG< ' Ts- Tc i ) /< Ts-Tco) ' > 

2040  A2-0tci/<  <  Ts-Tci  >"»t_GG<  <  Tj-Tci  )/(  Ts-Tco)  >  ) 

2050  A3-0tco/< ( Ts-Tco >»L0G< l Ts-Tci>/< Ts-Tco>  > > 

2060  Dlmtd'Lnitd»SQR(A1  '2  "A2  '2*A3  '2 ) 

2070  Dao  =  Uo»SGR( (Dqp/Qp)  2*< D imtd/Lmtd )  '2 > 

2060  M=Mu 

2090  Tl=< T*273.  i5>/273. !5 

2100  Ku-FNKuan 

2110  Ac-0.        !    INTERSCEPT   "ROM  SIEDER   PROGRAM 


135 


2120 

2i'JC 

2140 

2 '50 

2160 

2170 

2!  00! 

2190 

2200 

2210 

2220 

2240 

2260 

2270 

2290 

2290 

2300 

2310 

2320 

2330 

2  340 

2250 

2?&0 

2370 

2380 

2330 

24U0 

2^0 

2«20 

A4> 

2«30! 

2440 

2*50 

2460 

l  >  >) 

2470 

24P.0 

2*90 

2500 

2510 

2520 

2530 

";5^0 

2550 

2560 

2570 

2580 

2590 

2600 

2b  10 

2620 

2630 

2A40 

2650 

2660 

2670 

2680 

2t90 

2700 


LENGTH  OF  UNFTMNEf-  LEFT  PART  OF 


:ngth  OF  .'NFINNED  right  PART  OF 


UBE 

TUBE 


LI -.060325 

L2- .034925 

D1-. 01905 

D2-.  0:5675 

Pr"Cpu*Mw/Kw 

Muu*FNf1uu  <  T ) 

UNCERTAINTY  OF  INSIDE  HEAT-TRANSFER  COEFF. 

C*-1. 

Hi  =  <Ku/Di>»<Ci»R<?  ■  .8-Pr    .333»Cf-Ac) 

Dti»Q/(PI*Di»(L*Ll»Fel+L2»Fe2>«Hi) 

Cfc=<Mu.w/FNMuw<  F+Dt:  >>~.14 

IF    ABSMCfc-Cf  >/Cfc>>.01    THEN 

Cf-<Cf+Cfc>«.5 

GOTO   2200 

END    IF 

P1-PI*(Di+Dl > 

B1=(D1-Di)«PI*(Di+D: >-.5 

M1-<Hi-P1 /(kc-31 > >     .5 

P2-?I*<Di+D2> 

B2-(D2-Di)*PXDi+D2>-.5 

N2-<Hi*P2/iKc*B2) )    .5 

Fei  «FN7anh<M1  »L1  )/<M1  »»_1  > 

Fe2«'rN7anh<M2-L2>7CM2-t-2)  i 

Dtc*Q/fPI»Di»(L+L'  »Fe1  +  l2*-~<?2>»H>  > 

IF  A6S(<Dtc-Dti)/Dtc)>.01  ThEN  2200 

Dku- .  001 0  !  ERROR  IN  WATER  THERMAL  CONDUCTIVITY 

Dci-.0005  !  ERROR  IN  SIEDER-'ATE  COEFFICIENT 

Ddt».05  !  ERROR  IN  PRANDT.  NUMBER 

Dct-S.E-6 

A4  =  .  ia»Dcf  /Cf 

Dhi-hi»SQR<(Dku/Kw)  2*<Ddi/Di>  2+< .3»Dre/Re>  2+< . 333-Dpr/Pr )  2+<Dci/Ci >  2* 

UNCERTAINTY  OF  OUTSIDE  HEAT-TRANSFER  CHEF'. 
Rw  =  Do-lOG<Do/Di>/<2»>.c>  !  WALL  RESISTANCE 
Ho*  '/<  (  !/Uo)-<Do*L/(Di*<L+Ll-Fe1  +  L2-'re2>-Hi  >  )-Ru) 

Drw-Ru»SQR<  (Ddo/Do> '  2+  (  Dk  :c/Kc  )  "2  +  <  Ddo/  (  Do«L0G<  Do/Di  )  )  )  2*<Ddi /< Di «lOG(Do/D 
2) 
A5«1/Uo-Rui-<Do/<Di«Hi>> 
A6»Du.o/<Uo'  2»A5> 
A7-Dru/A5 

A8  =  <  <Do/(Di-Hi ) )»(Dhi/Hi ) >/A5 
PRINT 

Dno=Ho~>QR< A6  2+A7~2+A8 '2> 
CALCULATE  THE  7.  UNCERTAINTY  IN  Ho 
Prho-Oho- lOO/Ho 

CALCULATE  THE  7.  UNCERTAINTY  IN  REYNOLDS  NUMBER 
Prre-Dre»100/Re 

CALCULATE  THE  7.  UNCERTAINTY  IN  MASS  FLOW  RATE 
Prmf-Dm*"»lOO/Kf 

CALCULATE  THE  '!  UNCERTAINTY  IN  HEAT  TRANSFER 
Prqp=Dap*100/Qp 

CALCULATE  THE  '■'.    UNCERTAINTY  'IN  LMTD 
Prlmtd«Dimtd*100/'Lmtd 
CALCULATE  THE  >.   UNCERTAINTY  IN  Ru 
Prru=Dru»100/Ru 

CALCULATE  THE  7.  UNCERTAINTY  IN  OVERALL  HEAT  TRANSFER  COEF. 
Prac«Dao-100/Uo 

CALCULATE  THE  7.  UNCERTAINTY  IN  INSIDE  HEAT  TRANSFER  COEFF. 
Prru-Dhi*100/H: 
PRINT 
PRINT  USING  "10X. ••••UNCERTAINTY  ANALYSIS: 
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2710  PRINT 

2711  PRINT  USING  "10X.""    VARIABLE  PERCENT  UNCERTAINTY 

2712  PRINT 

2720  PRINT  USING  "1  OX  .""^ass  Flow  Rate.  Md  ,'"  .  ?D  .  2D  .  "  :?rmf 

2730  PRINT  USING  "1  OX  .  ""Reynolds  Ntunoer  .  Re  ""  .  3D  .  2D  .  "  :Prre 

2740  PRINT  USING  " 1  OX . ""Heat  Flux,  q  "" . 2 . 2D . " : Prop 

2750  PRINT  USING  "1  OX .""Log-Mean- Tern  Diff.  LMTD  ""  .  Z  .2D"  :Pr imtd 

2780  °RINT  USING  "1 0X.""Wal 1  Resistance.  Ru  ""  .  DD . 2D  .  "  :  Pr ru 

2770  PRINT  USING  "". OX  . ""Overai  1  H.T.C..  Uo  ""  .  DD  .2D  ,  "  :Pruo 

2780  PRINT  USING  "!0X. ""Water-Side  H.T.C..  Hi  "" . 2D . 2D . " :Prhi 

2790  PRINT  USING  "10X, ""Steam-Side  H.T.C..  Ho  "" .3D . 2D . "  ;Prno 

2300  END 

2810  DEF  FNMuu(T) 

2920  A-247.8/<T+133.1S> 

2830  Muu=2.4E-5*10~A 

2840  RETURN  Mum 

2850  FNF.ND 

28S0  DE~  -MTann(X> 

2870  P-EXP(X) 

2880  Q-EXP<-X> 

2890  Tanh-<P-0)/(P+ij) 

2900  RETURN  Tann 

2910  FNEND 

2920  DEF  FNKu<T! ) 

2930  Ku--. 32247+ Tt-<2.3395-T'«< I .8007-T I •< .5257  7- . 0  7344-T 1 ) >> 

2940  RETURN  Ku 

2S50  FMEND 

2980  DEF  "Nnu'(T> 

29T0  A-2*7.3/<TM33.15> 

2980  Mu-2.4F-5-l0  A 

2990  RETURN  Mu 

3000  FNEND 

30  10  DEF  FNRho' T) 

2020  Rho= 999. 52948+ T-<  .  01289-1 *< 5 . 4825 !3E-3-T» 1 .234147E-5) ) 

3030  RETURN  Rno 

3u40  FNF.ND 

3050  DEF  FNCdw<T> 

3060  Cpw-<4.2: i20858-T»<2.26326E-?-T«<4.4226lE-5  +  2.7'429E-7»T) ) ) - 1 000 

30  70  RETURN  Cpw 

3080  FNEND 

3090  DE^  FNTvsv<Emf  > 

3  UI0  COM  /Cc/  C(7> 

3110   T-C(C) 

3120  FOR  I=i  tq  ? 

3130  T-T+C(I>»Emf  I 

J!  40  NEXT  I 

3150  RETURN  T 

3160  FNEND 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS'. 


rile  Name : 
P r  e s  5 ur  e  Co n d  i  t  i o n  : 
Stean  Temperature 
hater  F  iou  Rate  ('/< ) 
Nater  Velocity 
heat  F  i  ulx 


f  *>  O  U  ?  c 

Vacuum    <  !  !    --  F'a) 
AS  .  /|7    (Deg   C> 
20.00 
i . 1 G    (n/s) 

2.111E+05    (W/m'2) 


UNCERTAINTY  ANALYSIS: 
VARIABLE 


:'ERCENT  UNCERTAINTY 


Mass  Floy  Rate.  Md 
Reynolds  Number,  Re 
Heat  Fiu.x .  q 

t_  -'  '^    I  '  T  CLI  I     r  ^  ■■•!    :_'  A  I   T   .    I .  <  14.' 

Ua.ll  Resistance.  Ru 
Overall  H.T.C. ,  Uo 
Wat e  r -S i de  H  .  T  .  C  .  .  Hi 
S  t  e  am  -  S  i  de  H  .  T  .  C  .  .  Ho 


5 

,44 

.48 

*. 

.74 

0 

o? 

2 

iS7 

1 

.75 

A 

.46 

>  -> 

,c.3 

Ijorq  rr.R  THE  UNCERTAINTY  ANALYSIS": 


'-lie  Name : 
p  r  e  s  s ur  e  Condition: 
Steam  Temperature 
Water  Floy  Rate  '*/> 
l-^at  er  Ve  1  oc  l  ty 
Heat  Flux 


~?2M/'2F 

Vacuum  (  !  1  k  Da » 

48.  ?8  (Deg  C  ) 

8G  .  00 
A  .  -' '  i  n '  <•  ) 
3.308E+Q5  (W/m  2) 


UNCERTAINTY  ANALYSIS 


VARIABLE 


PERCENT  UNCERTAINTY 


Mass  Floy  Rate.  Ma 
Reynolds  Number.  Re 
Heat  F  i  ux .  q 
Log-Mean-Tem  Diff,  LMTD 
Hail  Resistance.  Ru 
0  /eral  1  H.T.C,  Uo 
Water-Side  H.T.C. .  Hi 
Steam-Side  H.T.C. ,  Ho 


44 
57 


43 

E7 
66 
49 


2.31 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS 


r  j  i  e  'lame  : 
Pressure  Condition: 
3  team  Temperature 
Water  Fiou  Rate  (*/) 
Water  Velocity 
Heat  Fiux 


Atnbspher ic  (101  kPa) 
3S.33  (Dea  C) 
=   20.00 

I  .  IS  («n/s) 
-  6.502E+C5  (W/ni*2) 


UNCERTAINTY  ANALYSIS: 

VARIABLE 

Mass  F low  Rate.  Md 
Reynolds  Number.  Re 
Heat  F lux .  q 


lean- 


'em  Dif*  .  LMTD 


kail  Resistance.  Ru 
Gverai 1  H.T.C,  Uo 
Water-Side  H.T.C.  .  Hi 
Stearo-Side  H . T . C  .  .  Ho 


PERCENT  UNCERTAINTY 


r 

r,  p 

j 

jU 

i 

74 

•'1 

;"l  "7 

V 

U  . 

0 

1— 

67 

1 

74 

4 

47 

->c 

7  ? 

A  -r0R  THE  UNCERTAINTY 


DIALYSIS: 


F • ie    Name : 
Pressure    Condition: 
S  t  e  an     i  e  m  r-  e  r  a  t  u  r  e 
water    Flow    Rare    ( 7 > 
Hater    Veicci cy 
Heat    Fiux 


t-ttnospneric    HQi    fcPa) 


3  2  A  f.  8 

a  *. 

h  • 

?9.33    (Deg   C) 
S  u  .  0  0 
4 .  yfj    (n/s) 
■    1  .  164F+06    <W/m' 2> 


UNCERTAINTY  ANALYSIS: 

VARIABLE 

Mass  Fiou  Rate,  Md 
Reynolds  Number,  Re 
Heat  F i  ux ,  q 
Log-Mean- Ten  Diff .  LMTD 
Wall  Resistance,  Ru 
Overall  H.T.C. .  Uo 
Water-Side  H.T.C. .  Hi 
S team-Side  H.T.C.,  Ho 


PERCENT  UNCERTAINTY 

i  .44 
!  ,r>9 
0.48 
0.14 
2.67 

1  .50 
2.80 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


rile  Name : 
Pressure  Cond 1 t ion : 
Steam  Temperature 
Water  Flou  Rate  (V«> 
Water  Veioc 1 ty 
Heat  Fiu.x 


F39V83 

Vacuum  M  !  kPa) 

-   48  50  (Deq  O 
=   20.00 

1  03  <  n  /  <= ) 
=  2.579E+05  (W/m'2) 


UNCERTAINTY  ANALYSIS: 
VARIABLE 

Mass  Flow  Rate,  rid 
Reynolds  Number,  Re 
heat  F  lux  .  q 
Log-Mean- Ton  Diff .  LttTD 
Wall  Resistance.  Ru 
Over a  1 1  H.T.C. .  Uo 
Water-Side  H.T.C. .  Hi 
Steam-b  i  de 


H.T.C. .  Ho 


DERCENT  UNCERTAINTY 

5.42 
I  .45 

1  .73 

0  .23 
4.39 

1  .75 
4.43 

j  ?  ot 


DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


<-  i  i  e  ;iane  : 
Pressure  Condition: 
Steam  Temperature 
Water  Flou  Rate  < '/  ) 
Water  Ve  J  oc  1  r  -/ 
Heat  Flux 


F3SV83 

Vacuum    I  !  !    :.'  -"'a) 
=      48.41    (fieg  O 
=      80.00 

3.S2    (n/si 
=    4.440E-tG5    (W/n'2> 


UNCERTAINTY  ANALYSIS: 

VARIABLE 

Mass  Flow  Rate,  Md 
Reynolds  Number,  Re 
Heat  Flux .  a 
Log-Mean- Tern  Diff ,  LMTD 
Wail  Resistance.  Ru 
Overall  H.T.C. .  Uo 
Water-Side  H.T.C. .  Hi 
St  earn- Side  H.T.C,  Ho 


PERCENT  UNCERTAINTY 

:  .43 
!  .55 
0.50 
0.48 
4.  39 
.6'=* 
1  .47 
I  .69 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS 


File  Name : 
Pressure  Condition: 
Steam  Temperature 

Water    Flow   Rate    (7.) 
kiater    ;^e  loc  1  ty 
Heat    F  1  llx 


F39A98 

^nospSer  ir  (  10  •  kPa) 
99.9!  (Deg  C> 
=  20.00 

!  .04  (m/c) 
-  8.131E+05  <w/m' 2) 


UNCERTAINTY  ANALYSIS 

VARIABLE 

Mass  Flou  Rate.  Md 
Reynolds  Number,  Re 

Heat  F  i  u.x  .  q 


Log 


■Mi 


Ln-Tc-rr,  Diff  .  LtfTD 


Wal  I  Resistance  .  Rw 
Overall  H.7.C.  .  Uo 
Water-Side  H.T.C..  Hi 
St earn- Side  H.T.C.  Ho 


PERCENT  UNCERTAINTY 
5.39 


r 

,43 

1 

,  72 

j 

,07 

4 

,39 

1 

T-1 

1 

.      <      4_ 

4 

.42 

.  b 

.67 

DATA  FOR  THE  UNCERTAINTY  AHALVSTS: 


rile  Name  t 
Pressure  Condition: 
Steam  "emperature 
Water  z  iou  Rate  (7. > 
^ater  ^s  ioc  .  ty 
Hear  F  j  llx 


F 39 A 98 

Atmospheric  <  ! 0 t  kPa) 

H9.73  (Deo  C> 

80  .00 
-,  ,  j  4  (  m  /  5  ) 

1  .  409E+0*:  <W/n'  2> 


UNCERTAINTY  ANALYSIS: 
VARIABLE 

Mass  Flow  Rate,  Md 
Reynolds  Numoer,  Re 
Heat  Flux ,  q 
Log-Mean-Ten  Diff  .  LfiTD 
kail  Resistance,  Rui 
Overall  H.T.C. ,  Uo 
Water-Side  H.T.C. ,  Hi 
St  earn- Side  H.T.C.  Ho 


PERCEMT  UNCERTAINTY 

1  .42 
1  .54 
0.48 
0.  !5 
4  .  39 
.50 
1  .46 
1  .69 
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